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DIGITAL CODES FOR DESIGN AND EVALUATION OF CONVECTIVELY COOLED 
ROCKET NOZZLE WITH APPLICATION TO NUCLEAR-TYPE ROCKET 
by John E. Rohde, Rudolph A. Duscha, and George Derderian 

Lewis Research Center 

SUMMARY 

The design of the coolant passages of a convectively cooled nozzle is dictated by heat- 
transfer and flow considerations that are compatible with the overall rocket system and 
acceptable stresses within the nozzle structure. The method presented for designing this 
type of nozzle utilizes both a design and an evaluation computer program used in conjunc- 
tion. These digital computer programs are coded in FORTRAN IV for use on an IBM 7094 
computer. 

The design program calculates coolant passage dimensions for a specified gas-side 
wall temperature distribution, and the evaluation program calculates a gas-side wall tem- 
perature distribution for specified coolant -passage dimensions. In addition, coolant con- 
ditions are calculated, and all heat -transfer and fluid flow results are utilized to calcu- 
late coolant-passage tangential and longitudinal stresses. 

The programs are set up to utilize a range of heat -transfer correlations, different 
tube surface roughness conditions, and a tube splice. The evaluation program has the 
additional capability to handle radial or flat-plate heat conduction across the tube crown 
and the application of a coating to the gas side of the coolant passage. 

Application of these programs is shown for the design and the evaluation of a nuclear- 
type rocket nozzle operating at chamber conditions of 4000° R and 530 pounds per square 
inch absolute with hydrogen as both the coolant and the propellant. The wide range of 
available heat -transfer correlations for both the hot-gas and coolant sides and the conduc- 
tion model used across the tube crown result in maximum gas-side wall temperatures 
ranging from 1550° to 2250° R. The majority of the calculated coolant-passage stresses 
are well into the plastic range for the material considered (347 stainless steel). 


INTRODUCTION 


Contemporary designs being developed for high-thrust propulsion systems employ 
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convergent-divergent nozzles. For liquid-chemical-fueled and nuclear rockets, the noz- 
zle is convectively cooled by the fuel. 

The coolant-passage design for this nozzle is dictated by heat -transfer requirements, 
flow considerations that are compatible with the overall rocket system, and acceptable 
stresses within the nozzle structure. The final design is the result of many iterations 
that involve the associated parameters of heat -transfer and flow considerations. The 
complexity of these calculations and the numerous iterations involved require the use of 
a high-speed digital computer. 

Computer programs to perform these calculations are in use on a limited basis 
throughout the rocket -nozzle-development industry and at government agencies. However, 
there is an apparent lack of programs available for general use or being documented in a 
form suitable for publication. 

This report describes two computer programs developed at the Lewis Research 
Center and discusses their application. The programs are written in FORTRAN IV for 
use on an IBM 7094 computer. Usage of the programs will accomplish the required cal- 
culations necessary for the final design of a convectively cooled rocket nozzle operating 
at steady-state conditions. One program is employed for the design and the other for the 
evaluation of this type of nozzle. The report is divided into a discussion of the analytical 
formulation of the two programs and an application to a given design that also includes 
determination of the effects of variations in available heat -transfer correlations and heat- 
conduction model. 

The design program utilizes the desired gas wall temperature profile as an input and 
calculates the coolant-passage dimensions required to match it, while the evaluation pro- 
gram utilizes fixed coolant-passage dimensions as an input and calculates the resulting 
temperature profile. The preceding calculations are performed by an interrelated series 
of heat -transfer and fluid-flow equations. The state of the coolant as it progresses 
through the passage is also calculated. The results of the preceding calculations are 
used at each specified axial increment to determine the magnitude of the tangential and 
longitudinal stresses in the coolant-passage structure. The formulation of the coolant- 
passage stress analysis is described in appendix C by Rene E. Chambellan. 

The model chosen for the application is a convectively cooled 347 stainless-steel 
conical nozzle with a throat diameter of 4. 3 inches, a contraction ratio of 12, and an ex- 
pansion ratio of 8. This model is considered to be a scaled version of a nuclear -rocket 
nozzle. Because nuclear rockets utilize hydrogen for attainment of high specific impulse, 
hydrogen was chosen as the propellant and the coolant with equal flow rates of 16. 7 pounds 
per second. Chamber conditions are 4000° R and 530 pounds per square inch absolute, 
and coolant inlet conditions are 56° R and 1150 pounds per square inch absolute. 

Coolant -passage configurations for convectively cooled nozzles have been rectangu- 
lar, circular, elliptical, and U- or D-shaped, as shown in figure 1. The coolant -passage 
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geometry of each of these programs has 
been restricted specifically to the use of 
flow areas for U- or D-shaped tubes. This 
choice of geometry was made because U- or 
D-shaped tubes are currently widely used 
in the construction of nuclear -rocket -engine 
nozzles. 

In addition to the final design presented, 
the use of the available heat -transfer corre- 
lations leads to a wide range of predicted 
values for coolant -passage wall tempera- 
tures. These results for a fixed design can 
be an aid to revising the design to allow a 
slight amount of conservatism or safety 
margin. A study of these effects was made 
based on the design achieved for the previ- 
ously mentioned conditions, and these re- 
sults are also presented. 

DESCRIPTION OF COMPUTER CODES 

Figure l. - Coolant-passage configurations for convectively 

cooled nozzles. Basic Analytical Procedure 

The design program and the evaluation 

program are similar in that the same heat -transfer and pressure -drop equations are 
used. The basic difference between the programs is that the design program utilizes 
specified gas-side wall temperatures and determines the required coolant -pas sage geom- 
etry, while the evaluation program uses a given coolant-passage geometry and calculates 
the gas- side wall temperature. 

Design procedure . - The design procedure is similar to that discussed in reference 1 
and is as follows: 

(1) Specify the input data 

(a) Chamber fluid conditions, hot-gas-expansion data, flow rates, passage-wall 

thickness, gas-side contour dimensions, coolant inlet conditions, number of 
tubes, and material and fluid properties 

(b) Gas-side wall temperature 

(2) Calculate the hot-gas heat flux into the nozzle-coolant-passage wall 

(3) Calculate the coolant-side wall temperature 
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Wall thickness 


Heat flow 


(4) Calculate coolant -passage dimensions required to match specified wall tempera- 

tures 

(5) Calculate the coolant state 

(6) Iterate items (4) and (5) until convergence is obtained on the coolant state within 

prescribed limits 

(7) Calculate the coolant-passage wall stresses 

Limitations in the form of geometry and flow conditions are imposed on this procedure. 

Figure 2 illustrates the geometry of 
a U-tube coolant passage. The approxi- 
mate flow area consists of a semicircle 
of radius R, . combined with a trape- 

t, 1 

zoid of a height equal to the total passage 
height minus R .. (All symbols are de- 
fined in appendix A. ) The area of the 
semicircle of radius Rx * is considered 

Ly 1 

as being the minimum allowable passage 
flow area. If the area required for heat- 
transfer purposes is smaller than this, 
the minimum semicircular area and 
shape are used. The resulting tempera- 
tures are calculated by an evaluation 
procedure within the design program. 

A coolant velocity limitation in the 
form of a maximum allowable Mach number is also set to preclude choked flow and ex- 
cessive pressure drops. When this velocity limit is reached, the required coolant flow 
area not to exceed this limit is determined, and the resulting temperatures are calcu- 
lated by the same evaluation procedure as mentioned previously. 

Another limitation related to the maximum allowable Mach number is a limitation on 
the coolant static pressure. At high coolant velocities, conditions could exist, from one 
increment along the coolant passage to the next, so that all the coolant pressure could be 
lost. The reason is that in the established sequence of calculations, pressure drop is 
calculated prior to the Mach number calculation. Therefore, before the limiting Mach 
number routine increases the passage flow area, pressure conditions could be obtained 
that would give unrealistic results. This possibility can be avoided by setting a minimum 
coolant static pressure as a percentage of inlet pressure. 

When this limit is reached, at any station for which passage dimensions are being 
determined, the prescribed gas-side wall temperature at that station is too low and must 
be increased. Consideration is first given to the possibility of the calculation being made 
in the region of the nozzle where wall temperatures are decreasing, usually beyond the 
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throat in the direction of coolant flow. The temperature is reset equal to a temperature 
200° R lower than that existing at the previous station. Calculations are then repeated. 

If this station were in a region of the nozzle where temperatures would be still increasing 
at subsequent stations, a minimum pressure could again be reached and a second adjust- 
ment of wall temperature would be required. This time the temperature is set equal to 
that of the previous station. If the minimum pressure is reached a third time, the pro- 
gram will stop and print out the results for each station up to the troublesome one. The 
previously mentioned temperature adjustment routines in most cases should yield results 
that are sufficient to allow the calculations to continue for the rest of the nozzle. How- 
ever, if two adjustments are not adequate, the program is stopped and the existing results 
are examined. This procedure is used rather than attempting to anticipate and adjust 
within the program for any further unrealistic input values. If two adjustments of wall 
temperature are adequate to avoid reaching the minimum pressure, the Mach-number - 
area-adjustment routine could still cause an additional increase in the temperature, as 
described previously. 

Typically, the nozzle design resulting from the design program is refined by exterior 
adjustments that are then used as input to the design or evaluation program. These re- 
finements result from a trade-off between heat transfer, pressure losses, material se- 
lection, stresses, and fabrication complexity. The predominant refinement of the final 
result from the design program is the smoothing of the coolant -passage heights. For 
fabrication purposes, the pressure shell axial profile desired is a smooth continuous 
curve along the nozzle. 

Evaluation procedure . - The evaluation procedure is similar to the design procedure 
and is as follows: 

(1) Specify the input data 

(a) Chamber conditions, hot-gas-expansion data, flow rates, passage-wall thick- 

nesses, gas-side contour dimensions, coolant inlet conditions, number of 
tubes, material and fluid properties 

(b) Coolant-passage dimensions 

(2) Estimate the hot-gas heat flux into the nozzle coolant -pas sage wall 

(3) Calculate the gas-side wall temperature 

(4) Calculate the coolant-side wall temperature 

(5) Calculate the coolant state 

(6) Calculate the heat flux being added to the coolant 

(7) Iterate items (2) to (6) until a balance exists between the two heat fluxes within 

convergence limits 

(8) Calculate the coolant-passage wall stresses 

A limitation in the form of a minimum coolant static pressure is imposed on this 
procedure to avoid complete loss of all the coolant pressure. When this limit is reached, 
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at any station, an informative message is printed and the evaluation program is restarted 
with the inlet coolant pressure increased by 50 pounds per square inch. This limitation 
will usually allow the complete evaluation of a nozzle in one run, which could require 
several restarts. 

Variations of desi gn and evaluation proced ures . - Each program has the provision 
for incorporating the use of tube splices. In regions where the nozzle contour diameter 
is much greater than that at the throat, the tube radius to thickness ratio could be so 
large from a stress consideration associated with tube buckling that a tube splice would 
be required. By doubling the number of tubes, this ratio would then be reduced by a fac- 
tor of 2, and the potential for tube buckling would be reduced. The programs only con- 
sider the use of a bifurcation as being practical from a fabrication viewpoint. 

For the nuclear-rocket nozzle, the back or shell side of the coolant passage, which 
is normally considered a non-heat-transfer surface, does transfer into the coolant the 
heat generated within the nozzle pressure shell as a result of nuclear heating. Both pro- 
grams have a provision to account for this heat as an addition to the total amount of heat 
the coolant absorbs. However, the resulting temperature rise of the shell is not con- 
sidered in the stress analysis. In addition, both programs have provision to include 
thermal radiation as an addition to the hot-gas-side heat flux utilized. The radiation 
heat -flux is specified in the form of an input constant at each nozzle station. 

Provision is made in the evaluation program for the utilization of a coating on the 
gas side of the coolant passage. The program accepts the coating in the form of either 
a thermal resistance or a specific coating and thickness of coating at each section along 
the nozzle. The use of a coating is presented in the section Use of Thermal Barrier 
Coating. 


Heat-Transfer Equations Utilized 

To accomplish the heat balance through the coolant-passage wall, equations are re- 
quired for the convective heat -transfer coefficients on both the hot-gas side and the cool- 
ant side. References 2 to 9 discuss the extensive amount of analytical and experimental 
work that has been conducted in the area of heat transfer, which can be utilized for con- 
vectively cooled rocket nozzles. The following correlations are a comprehensive collec- 
tion selected to encompass the range of correlations available. However, the coolant- 
side correlations have been restricted to the extent that only correlations obtained for 
hydrogen were considered. 

Hot -gas side . - Typical of the hot-gas-side equations is the Bartz equation of refer- 
ence 2 

Nu f = C g Re®‘ 8 Pr 8 * 4 (1) 
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where = 0. 026 and where the properties are evaluated at a film temperature that is 
the average of the gas-side wall and the gas-static temperatures. 

Reference 3 has shown that this type of correlation is valid if appropriate values of 
Cg are used with properties evaluated at a reference condition, as described subse- 
quently. The authors of reference 3 indicated that C is a function of nozzle area ratio 
and could vary from their experimental values for other nozzle contours. They also indi- 
cated an exponent of 0. 3 for the Prandtl number. 

The Nusselt equation utilized for the computer programs is 

Nu f = C g Re f ' 8pr f ' 3 ( 2 ) 

The C values are supplied as input to the program. Figure 3 shows three distinct 

o 



Figure 3. - Various distributions of hot -gas-side correlation constant as function of nozzle area ratio. 


curves of as a function of nozzle area ratio. One is from reference 3, another is 
from reference 4, and the third has been determined from unpublished heat -transfer data 
and will be referred to hereinafter as the design curve. The curve from reference 4 is 
based on equation (1) with properties evaluated at a film temperature defined as the aver- 
age of the gas-side wall and the gas effective temperature, which is discussed in the 
following paragraph. Both the curve from reference 3 and the design curve are based on 
equation (2) with fluid properties evaluated at a reference condition, which is also dis- 
cussed in the following paragraph. Values of should be selected by taking into con- 
sideration the experimental conditions for which these values were determined. 

Reference 5 discusses the evaluation of properties for high velocity gas flow through 
the use of a reference enthalpy. If any variation of the hot-gas specific heat is neglected, 
a reference temperature for fluid properties evaluation of the same form as equation 
(10-16) in reference 5 can be used as follows 
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( 3 ) 


T ref * °' 5 T w, g + °- 28 T g, s + °' 22 T eff 


where 


T eff T g, s + R f ^ T g, tot ‘ T g, 
and the turbulent flow recovery factor 


R. 


t- V s 


ref 


( 4 ) 


For the range of conditions considered in this report, gas properties have been evaluated 
at reference temperature as indicated. 

The value of effective gas temperature is considered the driving temperature pro- 
ducing the heat transfer across the gas film in the relation 


(i) o = V T eff- T w >g ) 


( 5 ) 


The resulting equation for the gas-side heat -transfer coefficient, which is used in the 
programs, is 


C k , 
g g>f 


D 


H, g 


w \ D H P f 


-i0. 8 


- 


LvW M f p bJ 



Included in the preceding equation is the film Reynolds number, which is defined as 


Re^ = 


jv \ P H p f 

A f.J ^fPb 


where 


w 

A U 

and where b refers to a static condition. The term Pf/p^ is used to correct for the 
value of the bulk density appearing in the w/A^ term. 
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Coolant side . - Again, a type of Nusselt relation is used to calculate the heat -transfer 
coefficient on the coolant side as follows 

Nu f = C^Re 8 ’ 8 Pr°- 4 C X C 2 (6) 

where can be varied according to axial location or held constant. Properties are 
evaluated at a coolant film temperature, which is the average of the coolant- side wall and 
the coolant total temperatures. 

The coefficients C ^ and C 2 are optional corrections or modifications. The correc- 
tion Cj is that determined in reference 6 for turbulent flow and supercritical hydrogen 
at pressures to about 800 pounds per square inch absolute. The term accounts for the ef- 
fects of the extreme variations of the fluid properties across the boundary layer. It is 
normally applied while using as a constant and equal to 0. 0208. The relation for 
determining Cj is 


v 

C« = 1.0 + 0.01457 — (7) 

"b 

The existence of a secondary flow in curved pipes influences the heat-transfer coeffi- 
cient and friction factor, as shown in references 7 and 8. The term C 2 is a correction 
for the effect of coolant passage curvature on the heat -transfer process. Although this 
correction was derived from pressure-drop measurements, its adaptation to heat transfer 
was made because of the analogy between fluid friction and heat transfer. The following 
relation is employed that does not consider the effect of angular position as shown in ref- 
erence 8: 


C 


2 ' 



( 8 ) 


The original experimental friction factor results of reference 7 were published with the 
Reynolds number as the bulk Reynolds number. To allow some degree of conservatism, 
until better heat-transfer results for curvature effects are published, the Reynolds num- 
ber used for heat transfer in the C 2 term is that based on film conditions. 

Reference 9 indicates that, for straight tubes and supercritical hydrogen at pressures 
of 1000 to 2500 pounds per square inch absolute and turbulent flow, the relation 

Nu f = 0.021 RejMpr 8 ’ 4 


(9) 



is valid. However, it should be noted that equation (9) was derived from data where the 
wall to bulk temperature ratios ranged from 1. 5 to 11. 0. 

The resulting equation for the coolant-side heat -transfer coefficient, which is used 
in the programs, is 



with the option available for using Cj and/or Cg- Note that the film Reynolds number 
for the coolant side is defined exactly the same as the hot-gas film Reynolds number. 

Heat -transfer model. - To accomplish a heat balance across the coolant passage 
wall, a working heat -transfer model with some simplifying assumptions is required. Fig- 
ure 2 has illustrated the geometry of a U-tube coolant passage. The scalloped gas-side 
wall configuration is sufficiently complex so that the flow conditions cannot be defined 
readily in the braze region between the tubes. Previous experience leads to the observa- 
tion that the top of the tube crown has the highest temperature and that the wall tempera- 
ture diminishes at points closer to the tube braze. This temperature variation is further 
augmented by the fin effect of the tube walls, which conducts heat away from the tube 
crown. This leads to a complex two-dimensional heat -transfer problem in which the heat- 
transfer coefficients vary around the tube wall on both the inner and outer surfaces. The 
present state of development in the theory of nozzle design does not warrant a two- 
dimensional heat -flow analysis, because of its additional complexity. 

The objective of the heat-transfer model is to circumvent this complexity and yet ob- 
tain a realistic heat balance across the coolant passage wall. This balance is necessary 
for the determination of the metal wall temperatures and the coolant heat pickup. 

Meaningful metal wall temperatures can be obtained by considering the tube crown 
only, which is the region of maximum temperature. The tube portion exposed to the hot 
gas can be arbitrarily divided into two parts, according to the predominant mode of heat 
flow. One portion, the crown, transmits heat in a path essentially radial toward the tube 
centerline, whereas the remaining portion transmits heat in a skewed direction down 
toward the heat sink in the braze region where adjacent tubes meet. One -dimensional 
heat flow in the tube crown region can be considered radial or the area of heat flow for 
the outside and inside surface can be considered equal, as in the case of a flat plate. 

The true crown temperature is equal to or less than the maximum possible temperature 
based on pure radial heat flow, and the flat plate temperature is believed to set a lower 
limit for encompassing the true temperature. However, if experimental results estab- 
lish that temperatures well below the flat plate temperatures exist, a complex analysis 
as described previously could be justified. 
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The following equations are for one-dimensional radial flow of heat through a cylin- 


der 


T - T 
w, g w, £ 



where (q/A) Q is the heat flux entering the tube and is determined by 



~ V T eff 


- T. 


w,g 


) 


( 10 ) 


(ID 


The heat flux leaving the tube on the inside is 

= ('a') (12) 

A h [A 'o R t,i 

where the difference in the outside and inside areas is accounted for. 

The following equations are for the one -dimensional flow of heat through a flat plate 
of constant area 


w, g " w, £ 



(13) 


where (q/A) Q is calculated through the use of equation (11) and (q/A) i = (q/A) o< 

There is a significant difference between the temperature produced by the two meth- 
ods. As shown later in this report in the section Heat-conduction-model effect, the use 
of the radial heat -transfer equation results in a value of gas wall temperature 200° R 
greater than that obtained by the use of the flat-plate heat -transfer equation for typical 
conditions at the nozzle throat. Therefore, the use of the radial equation is more conser- 
vative and is used in both programs with the option of flat -plate heat flow being available 
in the evaluation program. 

Another area of concern, in completing the heat balance across the coolant passage 
wall, is the heat input to the coolant 
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( 14 ) 


n 


q = 



A., 
ht, n 

2-d, n 


where (q/A^ d n and A ht n are the outside heat flux and surface area over the incre- 
ment considered for a two-dimensional study. This heat flux would be the result of a two- 
dimensional analysis mentioned earlier as being too complex for use in the programs. 

The following equation can be used with the radial or flat-plate one -dimensional equa- 
tion as a basis 


q = 



A 


ht 


1-d, cr 


( 15 ) 


where the surface area used for heat transfer is 


A ht = M nR t, o, av e c, av 


where e is a heat -transfer-area correction. The value of e used does not ap- 

C f clV 

preciably affect the wall temperature, but it does affect the coolant heat pickup, as shown 
in the section Heat -transfer-area correction effect. A generally accepted value for 
e c av * s anc * va l ue I s used throughout this report, except where noted. 

The definition of the heat -transfer-area correction can be made in the following 
form: 


e 


c 



ft) 


2-d, n 


ht, n 


dL 


ht, sc 


cr 


If the two-dimensional study were performed at each station along the nozzle, a set of 

values for e could be established. By assuming that the established value of e at 
c i 

each station will remain constant for slight perturbations of conditions, the use of equa- 
tion (15) should yield results essentially equal to those of equation (14). This two- 
dimensional approach is recommended when attempting to match experimental coolant 
exit temperatures accurately. 


12 



Fluid Flow Equations Utilized 


Pressure drop . - Two factors, friction pressure drop and momentum pressure drop, 
comprise the complete coolant pressure drop calculated by the program. Any coolant 
passage inlet and exit losses are excluded and have to be considered for each separate 
design. 

The friction pressure drop is calculated from 


where G ay and D R> ay 
follows: 


ip Jr = 


2G av* 4t 


Ps, av^H, av^ 

are determined as average values over the increment AS as 


(16) 


G 


av 


w 2w 

%,av A iS, 1 + A f£,2 


^J^av = 4 ^ A f£, 1 + A f£, 2^ 
WP av ’ WP 1 + WP 2 


(17) 


(18) 


The value of p . is determined from an average static pressure and average static 

S, 3.V 

enthalpy over AS . . 

For smooth tube conditions, the friction factor is determined by equation (6 -8b) of 
reference 10, which is recommended for extrapolation to high values of Reynolds number. 
Equation (6-8b) is given as 


= 4. 0 log(Re av y/fj - 0. 40 

y/i 


(19) 


For rough tube conditions, the friction factor is determined by the equation 


VF 


-4. 0 log 



+ 



( 20 ) 


of reference 11, where e is the tube relative roughness. 
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These equations are used in the program, and the average Reynolds number is calcu- 
lated from 


Re 


av 


^av^H, av 




av 


( 21 ) 


where G ay , , and /u ay are average values over Af, evaluated as previously dis- 

cussed. 

A pressure drop curvature correction Cg from reference 7 is applied to the friction 
factor, similar to the heat -transfer curvature correction Cg that was applied to the heat- 
transfer coefficient. The equation 


C 


3 " 




( 22 ) 


is used where the Reynolds number is the same one as used for evaluating the friction 
factor. The term Cg is applied to the friction factor as follows: 

f c = C 3 f 

with f being used in equation (16). The range of validity for Cg and Cg is 

JhJ\ 2 a 


The momentum pressure drop is calculated from 


Ap 


mom 




(23) 


where ay is the average flow area over Af . 

Determination of coolant st agnation and static con ditions. - The inlet conditions, 
temperature and pressure, of the coolant are prescribed and are considered as being the 
stagnation or total temperature and the static pressure. For all practical purposes at the 
inlet to the coolant passage, total and static conditions are equal because of the low veloci- 
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ties usually existing there. These criteria have been adopted to facilitate the comparison 
of test data with analytical results. Generally, the inlet total temperature and static 
pressure would be measured because the instrumentation required is simpler. 

Initially, the static condition is assumed to be equal to the total condition, and itera- 
tions are required to arrive at the true conditions. From the relation 


Plot = p s + 


Pb v b 

2g 


(24) 


the total pressure at the inlet is calculated. The density used is evaluated at an average 
between the static and total conditions, which is assumed to be adequate for the low inlet 
velocities. The value of p^ and the prescribed value of T tot determine the total con- 
dition for obtaining the total enthalpy H^. By using this value for and the relation 


H s = H tot - 



(25) 


a value of static enthalpy H is obtained and when combined with the given pressure 

s 

yields the static condition. The iteration to arrive at an average density finally results 
in completely known static and total conditions. 

The total heat pickup by the coolant between the nozzle stations is based on an aver- 
age heat flux over the increment and heat -transfer surface area between stations: 


AH tot = 



(26) 


For the design of a nuclear rocket nozzle where heat is also generated within the nozzle 
pressure shell, provision is made for including this heat in the coolant enthalpy increase. 
The heat generation rate in Btu per cubic inch per second is specified as input data. The 
average value between stations is multiplied by the incremental shell volume, and it is 
assumed that all this heat enters the coolant. This value is then added to the heat trans- 
ferred to the coolant from the hot gas. Therefore, 

H to t 2 = H tot 1 + AH tot + Nuclear heat addition (27) 

is the enthalpy of the coolant at station 2. The static condition at station 2 is determined 
by using this enthalpy in equation (25) and the static pressure calculation described in the 
section Pressure drop. 
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The total condition of the coolant will be known when the iteration on total pressure 
is completed. For large differences between the total and static conditions, the varia- 
tions of the density term in equation (24) must be known. To avoid the problem of deter- 
mining the variation between those two conditions, a relation was derived for the varia- 
tion of total pressure, through the use of the energy equation and defining state equations 
The derivation is described in appendix B, and the result is as follows 


A Ptot = J Ptot, av AH tot ( 1 - - A Pfr| 


s, av , 


^tot, av T tot, av j 

p T ) 
x H s, av s, av / 


(B13) 


where all averages are arithmetic mean values from station 1 to station 2. It can be seen 
that, for the low velocity case, where static and total conditions are nearly equal, the 
change in total pressure is approximately equal to the friction pressure drop, which is to 
be expected. With this relation for the change in total pressure, the value of total pres- 
sure at station 2 is obtained, and the iteration is made by using the new values for total 
conditions in the Ap^ relation. Thus, the total fluid condition at station 2 is deter- 
mined. 


Coolant-Passage Stress Analysis Utilized 

The final design is based on considerations of heat transfer, pressure drop, and 
stresses developed in the structure. The temperatures of the coolant passage material 
alone are not necessarily indicative of safe operation; therefore, a stress analysis was 
incorporated into the programs. The stress calculations are made for each station, after 
the heat -transfer and pressure-drop routines have determined local temperatures and 
pressures. 

Since it is anticipated that thermal stresses due to the temperature gradients will be 
very high, the assumption that deformations will be elastic gives a qualitative character 
to the results. Because of the complexity of a plastic -elastic analysis, one was not in- 
cluded. Use of the stress-strain curves for the materials can give an insight into the de- 
gree of yielding that can occur. The calculated stresses can be of value for predicting 
the relative adequacy of a nozzle design. 

The idealized structure for the U-tube thick-shell configuration was assumed to con- 
sist of completely circular tubes supported in a bundle by a relatively thick shell, as 
shown in figure 4. Since the region of interest, the hot gas semicircular tube crown re- 
gion, is adequately represented by this configuration, it is assumed valid for the purpose 
intended in the programs. 
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Figure 4. - Tube and shell idealized model for stress calculations showing orientation of 
stresses. 


In general, the supporting conical shell is thick compared with the tube walls, which 
means that it will be much stiffer than the tubes. Therefore, the tubes will give negligi- 
ble resistance to deformation of the shell, and the tubes will deform as the shell dictates. 
No distinction is made for the fact that the flow passage diameter (which is the diameter 
of the tube crown in the U-tube configuration) varies along the nozzle length. Local con- 
ditions arising from bending of the shell at the throat and knuckle are neglected. The 
method used to analyze the nozzle stresses is discussed in appendix C, where the equa- 
tions used in the analysis are presented. 

The three pertinent calculated stresses shown in figure 4 are the total tangential 
stress at the inner and outer fibers of the tube crown and the average tube longitudinal 
stress. The tangential stresses consist of both bending and the membrane stresses at 
the crown. The longitudinal stress is caused by the net longitudinal strain in the tubes, 
due primarily to the large temperature difference within the material from the hot-gas 
side to the cold-shell side and the pressure forces acting on the shell. 

Another important structural consideration is the possible occurrence of buckling of 
the tube crown due to longitudinal stress. Since no adequate buckling analysis exists for 
this configuration, no attempt has been made to perform such an analysis for these pro- 
grams. ftowever, an analytical approach for similar situations is presented in refer- 
ence 12. This reference delineates that the longitudinal stress and the ratio of tube crown 
radius to tube thickness are important parameters with regard to tube buckling, and both 
should be maintained at as low a value as possible to prevent buckling. 
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The specific details of the use of both the design and the evaluation codes are de- 
scribed in appendix D. 


APPLICATION OF THE CODES TO A NUCLEAR-ROCKET-TYPE NOZZLE 

Design Procedure 

In the situation where a completely new nozzle design is required, various parameters 
would be known and fixed, a range could exist for some, and others would be open for ex- 
amination for effect. In the case of designing a nozzle for a nuclear rocket, the hot -gas 
conditions, that is, flow rate, and reactor exit temperature and pressure, would be set. 
For either flight or reactor test conditions, an exit area ratio would be established also. 

By knowing these values, an inside contour for the nozzle is determined for either a coni- 
cal or bell-shaped nozzle. The convergence region could be a subject for examination, 
such as determining the effect of convergence angle on the heat -transfer and stress re- 
sults. Usually there are governing factors such as overall length and the interface design 
between the reactor pressure vessel and nozzle pressure shell; therefore, it will be as- 
sumed that the complete gas-side contour and fluid condition would be set. 

On the coolant side, an exit temperature and pressure requirement would exist to 
allow a matching of system conditions. The coolant inlet temperature would be estab- 
lished according to supply conditions. The inlet pressure would be a variable up to a 
maximum value equal to the feed systems limitations. Therefore, two criteria associ- 
ated with the coolant to be achieved by the design would be a fixed minimum temperature 
rise and a pressure drop not exceeding the maximum allowable. The coolant flow rate 
for a nuclear nozzle would be essentially equal to the hot-gas flow rate. A slight devia- 
tion could exist to allow for a small percentage of bleed off from the main coolant supply 
to cool other parts of the system. 

The major parameter to be specified as input is the hot-gas-side wall temperature 
distribution. The actual distribution desired is difficult to specify until some initial re- 
sults have been obtained and the effect of these temperatures on the overall design deter- 
mined. As a start, a suitable approach is to specify a relatively low constant gas-side 
wall temperature and a Mach number limit of 0. 3. The Mach number, pressure, and 
minimum area limitations would then cause the program to recalculate new temperatures 
for the region where the low temperatures are unobtainable. A higher constant tempera- 
ture and a Mach number limit of 0. 3 could be specified and a new design obtained. The 
results of these two designs should reasonably bracket the final design desired. Typi- 
cally, the gas-side wall temperatures and related heat fluxes reach a peak near the throat 
and then drop to lower values in the regions fore and aft of the throat. 
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Nozzle conditions . - A scaled down configuration of a nuclear -rocket nozzle was 
chosen as the model for the example presented herein. Chamber conditions were set at 
4000° R and 530 pounds per square inch absolute, and coolant inlet conditions were set at 
56° R and 1150 pounds per square inch absolute. A throat diameter of 4. 3 inches was es- 
tablished. Hydrogen was the coolant and the propellant with equal flow rates of 16. 7 
pounds per second being used. Figure 5 shows the gas-side contour, axial station divi- 
sion of the nozzle and lists the basic dimensions. 



The coolant passage material selected was 347 stainless steel with properties as 
shown in figure 6. The initial selection of the number of tubes is made by selecting 
probable minimum values of tube wall thickness and tube crown radius that can be fabri- 
cated readily. A thickness of 0. 012 inch and a radius of 1/16 inch were assumed to be 
reasonable values based on existing fabrication information. By using these values, the 
throat diameter of 4. 3 inches, and a gap thickness for shim and braze material between 
tubes of 0. 010 inch, the number of tubes determined as an initial number for which to 
obtain results was 90. The use of other numbers of tubes is discussed later. 

The following assumptions were made to obtain the results for the design discussion 
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Figure 6. - Properties of 347 stainless steel as function of material temperature from curve fits of data. 

that follows. On the gas side, equation (2) was used in conjunction with the design 
curve of figure 3 (p. 7), with fluid properties being evaluated at reference temperature 
of equation (3). On the coolant side, because of the high coolant pressure range, equa- 
tion (9) was used. 

The program calculations were made at each of the 25 axial stations into which the 
nozzle was divided, as shown in figure 5. Because the calculations assume a linear rela- 
tion in each variable from one station to the next, the lengths between stations are not 
uniform. In regions where small changes in variables occur, relatively large distances 
were chosen. In the region near the throat where variables are rapidly changing, incre- 
ments were chosen small enough so that the assumption of linear changes would still be 
valid and not lead to erroneous results. 

Temperatures and heights of coolant passages . - A temperature of 1300° R and an- 
other of 1700° R were chosen as the initial constant values of gas -side wall temperatures 
for two design runs. Figure 7 shows the resulting temperature distributions produced by 
the program limitation features. The constant temperature of 1300° R was achieved up 
to station 8. At station 9, the Mach number limit of 0. 3 was exceeded, so that the pas- 
sage area was increased to just produce this flow condition. This resulted in a wall tem- 
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perature of 1514° R. This Mach number limitation was the governing factor at each suc- 
cesive station up to and including station 18. The required increase in flow area resulted 
in the calculated temperatures, as shown in figure 7(a). At station 19, the area required 
to produce 1300° R was smaller than the limiting semicircular minimum area. The use 
of the minimum area resulted in a temperature of 1809° R. This minimum area limita- 
tion was the governing factor at each successive station up to and including the final one 
and also resulted in temperatures higher than the specified input value of 1300° R. 

For the next computer rim, the constant temperature of 1700° R was maintained up 
to station 10. At station 11, the Mach number limit again prevailed and continued to pro- 
duce temperatures in excess of 1700° R up to station 18. For stations 19 to 21, the mini- 
mum area limitation governed and also produced temperatures greater than 1700° R. 
Conditions for stations 22 to 25 permitted the attainment of the specified 1700° R. 

The third temperature distribution, shown in figure 7(a), is that resulting from an 
input of 1700° R for the first 20 stations and 1300° R for the remaining five. Therefore, 
the results for the first 20 stations are identical to those for a constant input of 1700° R 
for the entire nozzle. At station 19, the minimum area limitation governed and prevailed 
for the rest of the remaining stations resulting in a temperature distribution similar to 
that for a constant input of 1300° R but higher by about 50° R at each station. 


TABLE I. - COOLANT TEMPERATURE RISE AND PRESSURE 
DROP FOR THREE TEMPERATURE DISTRIBUTIONS 


Total coolant 

Temperature input, °R 

conditions 

1300 

1700 

1700 and 1300 

Temperature rise, AT, °R 

94.3 

86.4 

87.8 

Pressure drop, Ap, psi 

180.3 

158.1 

160.5 


Table I gives the coolant temperature rise and pressure drop for each rim. The 
1300° R case produces the highest coolant temperature rise due to the cooler wall tem- 
peratures and resulting greater heat flux. This case also has the greatest pressure drop 
due predominantly to the smaller tube flow areas required to produce the cooler tempera- 
tures. The combined 1700° and 1300° R case with wall temperatures appreciably cooler 
than those of the 1700° R case at the last three stations produced no appreciable increase 
in coolant pressure drop over the 1700° R case. 

Figure 7(b) shows the axial distribution of total passage heights for each case. The 
combined 1700° and 1300° R case heights are the same as those of the 1700° R case up to 
station 21 and then are the same as those of the 1300° R case for the rest of the length. 
From the results indicated in figure 7 and table I, it is apparent that in the cylindrical 
part, upstream from the convergent region, the use of minimum areas at stations 21 to 25 
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(indicated by the constant value of total height) yields the most acceptable results. In the 
throat region from stations 8 to 18 either case results in irregular coolant passage 
heights. 

Starting with the results from the design procedure, the first step is the smoothing 
of the coolant passage heights from one station to another. These heights combined with 
the fixed gas-side wall diameters determine the inside shell contour. This contour must 
be selected to produce a smooth transition from one section to another without any dis- 
continuities and must also be based on matching the preliminary design coolant passage 
heights as closely as possible. The selected shell contour with resulting coolant passage 
heights is then checked through the use of the evaluation program, to determine if un- 



to) Total coolant-passage height. 

Figure 8. - Final results obtained through use of design and evaluation programs. Case 1. 
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acceptable material temperatures, coolant exit temperature, or coolant pressure losses 
are produced. This procedure is repeated until an acceptable shell contour is obtained. 

The final coolant passage conditions resulting from the aforementioned procedure 
are referred to as those of case 1 and are shown in figure 8. Figure 8(a) shows the final 
hot-gas wall temperature distribution and figure 8(b) shows the total passage height vari- 
ation along the nozzle. The coolant temperature rise is 87. 4° R and the pressure drop 
is 116.8 pounds per square inch. The gas-side wall temperature varies from approxi- 
mately 1500° R at both the divergent and convergent regions to 2250° R in the throat re- 
gion. The design Mach number limit of 0. 3 could be increased to allow smaller passage 
areas that would reduce this calculated temperature of 2250° R. However, this would 
also result in an increased pressure drop thereby possibly requiring an increase in cool- 
ant inlet pressure. 

This design was accepted as is, because conditions in the coolant passage of a nozzle 
are such that this should be conservative and actual temperatures should be lower. The 
correlations used on the coolant side did not include any curvature effect at the throat or 
any effect of coolant property variation due to the extremes existing between the wall 
temperature and the coolant temperature. Later it will be shown how the variations 
available produce a range of predicted temperatures. 

Coolant -pas sage stresses . - Figure 9 shows the stresses for these cases. It can be 
seen that basically they follow the contour and temperature distribution profiles. In all 
cases the maximum values for each specific stress are nearly identical. By referring to 
the yield strength curve of figure 6 (p. 20), it is apparent that almost all resultant calcu- 
lated stresses are in the plastic region. Therefore, conventional engineering criteria, 
based on yield strength, for predicting failure cannot be used. As the design becomes 
more nearly fixed, the conditions are severe enough to demand a detailed analysis of what 
really can be expected from a structural aspect. The predominant factor causing the high 
values of stress is the large difference in temperature from the hot-gas side to the cold 
shell. For the conditions being considered, it is difficult to reduce the large stresses 
through geometry modifications alone. However, if a thermal -barrier type of coating 
were put on a given design, this difference in temperature would be reduced markedly 
with a corresponding decrease in stresses and coolant exit temperature. The results of 
putting a coating on this design is discussed in the section Use of Thermal Barrier Coating. 

Variation in number of coolant tubes. - When the final temperature distribution of 
figure 8(a) was achieved, it was used for two more sets of design calculations, one for 
80 tubes and the other for 100 tubes. The 100-tube case results in comparable values for 
temperature rise and pressure drop of 86. 7° R and 126. 6 pounds per square inch, re- 
spectively. Although the gas-side wall temperature distribution achieved was identical 
to the 90-tube case, the minimum tube radius at the throat of 0. 0525 inch appears to be 
approaching a fabrication limitation. The 80-tube case results in higher temperatures in 
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{a) Inside tangential stress. 
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(c) Longitudinal stress. 

Figure 9. - Coolant -pas sage stresses as influenced by gas-side wall temperature input for 
design program. 

TABLE n. - GAS -SIDE WALL TEM- 


PERATURE FOR DIFFERENT 
NUMBER OF TUBES 


Station 

Temperature, °R 

90 Tubes 

80 Tubes 

19 

2000 

2005 

20 

1800 

1923 

21 

1734 

1885 

22 

1692 

1846 

23 

1638 

1795 

24 

1578 

1735 

25 

1509 

1664 


the cylindrical region because the minimum area re- 
striction for the wider tubes results in larger flow 
areas than those of the 90 -tube case. Table n shows 
the comparison of these temperatures for stations 19 
to 25. Only the 80-tube- and 90-tube-case tempera- 
tures are shown, because as noted previously, those 
for the 100 -tube case were identical to those of the 
90-tube case. The higher temperatures of the 80-tube 
case are not desired. Therefore, of the three different 
numbers of tubes, the value of 90 appears to be an ade- 
quate choice. 
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Effect of Assumptions on the Design 


This section illustrates how the evaluation program is used to determine the possible 
variations in predicted conditions during steady-state operation for a fixed design. These 
variations are produced by varying the correlations used and assumptions made during 
the design phase. If any of these variations produce an excessive wail temperature, pres- 
sure drop, or temperature rise, the results must be analyzed with the possibility of re- 
design or the addition of a coating. If the decision is made to redesign, the nozzle is sent 
back to the design program with the assumptions that proved to be troublesome being in- 
corporated into the design for conservatism. 

In the sections Heat-Transfer Equations Utilized and Fluid-Flow Equations Utilized, 
the following variations were presented: 



(b) Coolant correlation, Nu { - 0.0208 Re f °- 8 Pr f °- \l . 0 + 0.01457 v w /u b ). 
Figure 10. - Effect of gas-side correlation on gas-side wall temperatures. 
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<d) Inside surface heat flux. 

Figure 11. - Typical variations produced by use of different heat- 
transfer correlations on coolant side. 


(1) Three curves of values 

with two hot-gas heat -transfer 
correlations 

(2) Two coolant heat-transfer cor- 

relations 

(3) A curvature correction on the 

coolant-side heat-transfer 
correlation Cg 

(4) Two modes of heat conduction 

through the tube crown 

(5) A heat-transfer-area correc- 

tion for the heat added to the 
coolant 

(6) Two types of wall -surface con- 

ditions for friction pressure 
losses 

(7) A curvature correction on the 

friction factor Cg 

Heat -transfer -correlation effects . - 
Figures 10 and 11 show the effect pro- 
duced by the variations of items (1) 
and (2). The case numbers in the fig- 
ure are listed in table III with some 
of the significant resulting differences. 

Figures 10(a) and (b) show the ef- 
fect on the gas-side wall temperatures 
produced by varying the hot-gas-side 
correlation while using identical 
coolant-side correlations, respectively, 
as indicated in the figures. The hot- 
gas-side correlation and/or Cg distri- 
bution that is selected can be seen to 
have an appreciable effect on the wall 
temperature. An overall deviation of 
about 200° R exists between the maxi- 
mum temperatures of the three curves 
of each figure. 

Figure 11(a) is a replot of the tem- 
perature distributions of cases 1 and 5. 
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TABLE m. - HEAT-TRANSFER CORRELATIONS AND RESULTS FOR CASES EVALUATED 



Design 
Ref. 3 
Ref. 4 
Ref. 3 
Design 
Ref. 4 


Gas film temperature 


T f = T ref 
T f =T ref 

T f = ( T eff + T w,g )/2 
T f = T ref 
T f = T ref 

T f = ( T eff + T w,g )/2 


Coolant equation 


Nu f = 0.021 Re®' 8 Pr®‘ 4 
Nu f = 0.021 Re®’ 8 Pr®- 4 
Nu f = 0.021 Re®’ 8 Pr®' 4 
Nu f = 0. 0208 Re®' 8 Pr®’ 4 Cj 
Nu f = 0. 0208 Re®’ 8 Pr®‘ 4 Cj 
Nu f = 0. 0208 Re®' 8 Pr®’ 4 Cj 


Exit total 
tempera- 
ture, 
°R 

143.4 

137.3 

150.6 

155.7 
166.9 

181.3 


Exit total 
pressure, ] 
psia 

1043 
1050 
1037 
1029 
1012 
993 


see appendix D for program revisions necessary to run cases 3, 6, and 7. 


This figure is a comparison of the effect of identical hot-gas-side correlations and indi- 
cates that the effect of the different coolant-side correlations on wall temperature for 
these conditions is more significant. The deviation between maximum temperatures is 
about 500° R and is also the same for the other two pairs of curves shown in figures 
10(a) and (b). 

Figures 11(b) and (c) show a comparison of the coolant total pressure and tempera- 
ture distribution, respectively, for the same two cases, 1 and 5. These figures and ta- 
ble m indicate that the combination of hot-gas-side and coolant-side correlations has a 
significant effect on coolant conditions. For example, between the combination of cases 
2 and 6, the coolant exit temperature varies from 137° to 181° R and the coolant exit 
pressure varies from 1050 to 993 pounds per square inch absolute. 

The resulting variation of inside surface heat fluxes for cases 1 and 5 is shown in 
figure 11(d). The maximum heat flux varies from 14. 0 to 18.8 Btu per square inch per 
second in the throat region and drops to approximately 4 and 2 Btu per square inch per 
second in the convergent and divergent sections, respectively. Although not shown, the 
outside surface heat fluxes vary considerably from the inside surface heat fluxes for the 
case of radial heat flow. For example, the maximum heat fluxes near the throat vary 
from 11. 7 to 14. 0 Btu per square inch per second and 15. 7 to 18. 8 Btu per square inch 
per second for the outside and inside heat fluxes of cases 1 and 5, respectively. 

These variations indicate that the use of the design assumptions of case 1, design 
Cg curve, and the coolant correlation of reference 9, produce high wall temperatures 
on the average but relatively low coolant heat pickup. The remaining comparisons of 
items (3) to (7) are made by using the design assumptions of case 1 as the standard 
case. 
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Figure 12. - Gas-side wall temperature variations produced by curvature effect on heat 
transfer. Case 1. 


Curvature heat -transfer -correction effect . - Figure 12 shows the variations pro- 
duced in the wall temperature by the Cg correction in the localized regions of axial gas 
wall curvature. The heat -transfer enhancement due to curvature in the throat region has 
reduced the wall temperatures about 250° R. However, as mentioned previously the angu- 
lar effect is not considered in this correction, and that is why this correction does not 
produce a smooth gas wall temperature distribution. 

Heat-conduction-model effect . - Figure 13 shows the differences that exist if the 
heat -transfer calculations are based on either one -dimensional radial heat flow across 
the tube crown utilizing cylindrical geometry or normal heat flow utilizing flat-plate 
geometry. Figure 13(a) shows a reduction of approximately 200° R in maximum gas- 
side wall temperatures due to the use of flat-plate heat -conduction equations across the 
crown of the tube. Figures 13(b) and (c) show the resulting differences in coolant total 
pressure and temperature as being relatively small. 

Heat -trans fer -area correction effect . - Figure 14(a) shows the effect on the coolant 
heat pickup produced by a variation from 0. 75 to 0. 90 in the area correction This 
variation produces a significant increase in the coolant exit temperature of 21° R while 
producing only a 125° R variation in wall temperatures in the chamber region, as shown 
in figure 14(b). 

Coolant tube roughness and curvature effects . - Figure 15 shows the singular and 
combined effects of roughness and axial gas wall curvature on the coolant pressure. For 
the nozzle under analysis the rough tubes produced a decrease in exit pressure of 
53 pounds per square inch and the curvature reduced the exit pressure by 18 pounds per 
square inch. The value of tube relative roughness e used was 0. 00006 inch, the value for 
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Figure 15. - Effects produced on total pressure by tube roughness and curvature. Case 1. 


commercially drawn tubes. 

Combined effects . - As an additional check, the nozzle design is evaluated with the 
assumptions that will produce the maximum coolant heat pickup and pressure loss in the 
coolant passage, the lowest possible wall temperatures reasoned to exist in the nozzle, 
and the anticipated range of coolant passage stress values. 

Coolant: To achieve the maximum coolant heat pickup and pressure loss in the cool- 
ant passage the following assumptions were made: 

(1) Flat -plate heat conduction applies 

(2) A fluid-properties variation correction on the convective heat -transfer coefficient 

Cj applies 
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(3) The maximum Cg curve applies 

(4) The tubes are rough 

(5) A curvature correction on the friction factor Cg applies 
These assumptions will be referred to as case 7. 

Figure 16(a) shows the significant drop in wall temperatures by comparing the result- 
ing gas-side wall temperatures, case 7, against the designed values, case 1. As can be 
seen from figures 16(b) and (c), respectively, the bulk coolant temperature has increased 
by 41° R and the exit pressure has decreased by 210 pounds per square inch. However, 
to see the most significant effect, comparisons of the coolant velocity in each case must 
be made as shown in figure 17. This comparison gives an indication that the coolant flow 
may be in danger of choking. The maximum velocity has increased from 760 to 1100 feet 
per second and the maximum Mach number has increased from 0. 3 to 0. 5. In addition, 
this trend can be further augmented by increasing the area correction above 0. 8. 

Also, shown in figure 16 are the curves for minimum heat pickup and pressure drop, 
case 2. This case results in relatively negligible deviations of coolant conditions from 
the design, case 1, but produces an anticipated range of coolant conditions when combined 
with the maximum heat pickup and pressure drop curve, case 7. 

Wall temperatures: To achieve the minimum wall temperatures thought to exist in 
the nozzle, the following assumptions were made: 
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Figure 18. - Comparison of design gas-side wall temperatures with expected minimum 
gas-side wall temperatures. 


(1) Flat-plate heat conduction applies 

(2) A fluid-properties variation correction on the convective heat -transfer coefficient 

C x applies 

(3) The minimum curve applies 
These assumptions will be referred to as case 8. 

The resulting wall temperatures of case 8 are shown in figure 18 and, when com- 
bined with the temperatures of case 1, produce a range of anticipated wall temperatures. 
This range includes a variation from 1550° to 2250° R in the predicted maximum gas wall 
temperature. However, this range does not include the effects of coolant passage curva- 
ture on heat transfer that would further reduce the maximum wall temperature. 

Coolant passage stresses: As the coolant passage stresses are primarily a function 
of wall temperatures, the stress comparison will be made by using the maximum wall 
temperature range of cases 1 and 8. Figure 19 shows the anticipated range of stresses 
produced by the uncertainty in the heat-transfer analysis. It should be noted that a con- 
siderable drop in the magnitude of the stresses is produced by the lower material tem- 
perature of case 8. 
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(c) Longitudinal stress. 

Figure 19. - Variations in tube stresses produced by heat -transfer correlations. 



Use of Thermal- Barrier Coating 


The evaluation program has the provision for the addition of a coating to the gas side 
of the coolant passages. This coating would be added to an existing nozzle for uprating it 
or to a design with excessively high wall temperatures. Figures 20 and 21 show the ef- 
fects of two different thermal resistances on the designed nozzle. Figure 20 shows the 
marked decrease in the metal wall temperature and the marked increase in the gas-side 
coating wall temperatures. Figure 21 shows the decrease in all the tube wall stresses 
produced by the thermal resistance. However, the stress analysis presented in this re- 
port does not consider the coating. Although the tube stresses are reduced, the stresses 
in the coating and mismatch due to relative thermal expansion between tube and coating 
may still be excessive and would require further analysis. 

Figures 19 and 21 show that the variation in the predicted wall temperatures pro- 
duced a greater reduction in the stresses than the thermal resistance. However, it may 



Figure 20. - Gas wall and gas-side metal wall temperature profiles produced by addition of 
thermal barrier to designed nozzle for case 1. 
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(c) Longitudinal stress. 

Figure 21. - Tube stress variation produced by addition of thermal barrier to designed 
nozzle for case 1. 


be observed that the temperature reduction caused by the thermal resistance would also 
be superimposed on the reduction due to the use of the assumptions that produced varia- 
tions in the predicted wall temperatures. 

Off-Design Conditions 

The evaluation program can also be used to determine the conditions that exist at 
off -design points and the maximum extent to which a nozzle could be uprated based on 
some given criterion of failure. Figure 22(a) shows the high wall temperatures that re- 
sult when the nozzle chamber pressure is increased to 1000 pounds per square inch 
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(a) Gas-side wall temperature profiles. 



Station 


(b) Inside surface heat flux. 

Figure 22. - Results for uprated conditions of designed nozzle. Chamber pressure, 
1000 pounds per square inch absolute. 



absolute from the design value of 530 pounds per square inch absolute and the coolant in- 
let pressure is increased to 1500 from 1150 pounds per square inch absolute. The spread 
of 500° R between the wall temperatures of the two uprated cases indicates potential fail- 
ure of the nozzle if the correlations of case 1 were to apply. Figure 22(b) shows the range 
of inside surface heat fluxes that result from the uprating of conditions. 


General Remarks 

The previous discussion has been written in a general manner to give an insight into 
the approach to be used in obtaining a nozzle design. Each specific design has its own 
requirements that have to be obtained through a complete analysis whereby many parame- 
ters such as tube number, tube thickness, tube material, coolant inlet pressure, and 
possibly coolant inlet temperature and flow rate would be varied. For the purposes of 
this report, a detailed discussion describing all of the permutations possible leading to 
a final design was not deeme&hecessary. 


CONCLUDING REMARKS 

The large number of perturbations leading to any given final design for a nozzle 
coolant passage substantiates the requirement for the use of a design and evaluation pro- 
gram in conjunction, as shown in this report. Each program complements the other and 
together they produce the design, initial evaluation, and final resulting evaluation of fea- 
sible convectively cooled U- or D-tube nozzles with hydrogen as both the coolant and the 
propellant. 

The programs are set up to allow the option as to the heat -transfer correlations uti- 
lized, the condition of tube roughness employed, and the use of a tube splice. Additional 
options of radial or flat plate heat conduction across the tube crown and adding a coating 
to the gas side of the coolant passage exist in the evaluation program. 

The design program is used to determine a feasible nozzle coolant passage by using 
certain basic assumptions, fabrication limits, and a physical limit. The fabrication lim- 
its used were wall thickness of 0.012 inch and a minimum inside coolant tube radius of 
1/16 inch. The physical limit used was a Mach number limit of 0. 3 on the coolant flow. 
These limits were combined with reasonable wall temperatures and minimum coolant 
pressure loss requirements to produce a realistic design. 

The results presented are for a U-tube nozzle with a throat diameter of 4. 3 inches, 
a contraction ratio of 12, an expansion ratio of 8, and a total length of 31 inches. The 
nozzle is to operate at a chamber condition of 4000° R and 530 pounds per square inch 
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absolute and coolant inlet condition of 56° R and 1150 pounds per square inch absolute 
with hydrogen as both the coolant and the propellant. 

The final contour determined by the design program for the nozzle coolant passage 
was smoothed, and the effects of varying the basic assumptions were then checked by the 
use of the evaluation program. These variations include the following items: 

(1) Three curves of values with two hot-gas heat -transfer correlations 

( 2 ) Two coolant heat -transfer correlations 

(3) A curvature correction on the coolant-side heat-transfer correlation 

(4) Two modes of heat conduction through the tube crown 

(5) A heat -transfer-area correction for the heat added to the coolant 

(6) Two types of wall-surface conditions for friction pressure losses 

(7) A curvature correction on the friction factor 

The variations of the hot-gas and coolant -heat -transfer correlations employed pro- 
duced an overall predicted maximum gas-side wall temperature variation of 200° and 
500 ° R, respectively. The two modes of heat conduction across the coolant tube crown, 
radial and flat plate conduction, produced a predicted maximum gas-side wall tempera- 
ture variation of 200 ° R. Gas wall axial curvature produced localized regions near the 
throat with 250 ° R variations in the predicted gas wall temperatures. The combined ef- 
fects produced a variation in the predicted maximum gas-side wall temperature from 
1550 ° to 2250 ° R. 

The variations in the heat -transfer correlations also produced a significant variation 
of 44° R in the coolant exit temperature and 57 pounds per square inch in the exit pres- 
sure. The effect of using a tube roughness equal to commercially drawn tubes resulted 
in an increase in pressure drop of 53 pounds per square inch, over the original pressure 
drop of 117 pounds per square inch for smooth tubes. 

In general, calculated tangential and longitudinal stresses for the coolant passage 
wall are well into the plastic region. Therefore, their usefulness for predicting failure 
is doubtful, but they can be used to evaluate various designs on a relative basis. 

In addition, the scope of the final evaluation can include the following: variations in 
operating conditions such as uprating the nozzle; use of coatings; analyzing nozzle fail- 
ures; and matching experimental nozzle test data to determine the validity and applicabil- 
ity of correlations and assumptions. Examples of the uprating and coating addition were 
presented for the chosen application. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 8, 1966, 

122 - 29 - 07 - 07 - 22 . 
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APPENDIX A 


SYMBOLS 


A 

area 

J 

mechanical equivalent of heat 

A* 

nozzle -throat area 

k 

thermal conductivity 

C 

numerical coefficient 

L 

cone slant height from previous 

C 1 

convective heat -transfer coef- 


station 

ficient correction for fluid- 
properties variation 

AH 

linear distance along coolant- 
passage wall between stations 

C 2 

convective heat -transfer - 
coefficient correction for 

M 

bending moment acting on tube 
crown 


tube curvature 

N 

membrane force per unit length 

^3 

friction -factor correction for 


in circumferential direction 


tube curvature 

NS 

membrane force per unit length 

c p 

specific heat at constant pres- 


in longitudinal direction 

D 

E 

e 

f 

sure 

diameter 
Young’s modulus 
relative roughness of surface 
friction factor for straight 

N t 

Nu 

Pr 

P 

number of coolant tubes 

Nusselt number 
Prandtl number 
pressure 

f 

tubes 

friction factor for curved 

APfr 

friction static-pressure drop 
between stations 

c 

G 

tubes 

mass flow per unit cross- 
sectional area 

^mom 

q 

momentum static -pressure 
drop between stations 

local heat flow rate 

g 

gravitational conversion 

R 

radius 


factor 

R c 

radius of curvature 

H 

enthalpy 

R f 

recovery factor 

AH 

enthalpy increase between 
nozzle stations 

Re 

Reynolds number 

h 

convective heat -transfer coef- 
ficient 

S 

cone slant height from present 
station 
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s 

entropy 

fi 

flow 

T 

temperature 

g 

gas 

t 

thickness 

H 

hydraulic 

V 

velocity 

h 

hoop 

W 

deflection force acting on tube 

ht 

heat transfer 

WP 

wetted perimeter 

i 

inside 

w 

mass weight of flow 

int 

internal 

a 

coefficient of thermal expansion 

SL 

liquid or coolant 

6 

tube deflection 

long 

longitudinal 

e 

strain 

M 

meridional 

e c 

heat -transfer surface-area cor- 

m 

mean 

M 

rection factor 
dynamic viscosity 

max 

met 

maximum 

average between inside and outside 

V 

"n 

kinematic viscosity 
Poisson’s ratio 

n 

metal-wall temperatures 
division of cross section into n 

r' 

p density 

a stress 

cp angle between nozzle side wall 

and nozzle axis 

Subscripts: 

0 

P 

ref 

nodes for two-dimensional heat- 

th 

conduction study; also n n node 
outside 
pressure 
reference 

A 

axial 

s 

static 

av 

b 

average 

bulk 

sc 

sh 

semicircular 

shell 

bend 

bending 

T 

temperature 

com 

combined 

t 

tube 

cr 

tube crown 

tan 

tangential 

d 

deflection 

tot 

total or stagnation 

e 

external 

w 

wall 

eff 

effective 

y 

yield 

f 

film 

l 

previous station 
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2-d two dimensional 


2 present station 
1-d one dimensional 9 circumferential 



APPENDIX B 


STAGNATION PRESSURE STATE 


For a constant density fluid, the following relation for stagnation pressure is ade- 
quate: 


Ptot = Ps 



(Bl) 


However, when the density varies with both temperature and pressure, the evaluation of 
the proper density value to use becomes involved. To circumvent this problem, a differ- 
ent approach to stagnation pressure determination was used. 

The variation of stagnation pressure in a fluid process can be obtained from the en- 
ergy equation. The process of attaining the stagnation state is isentropic, and, there- 
fore, 


From the definition of 


ds, , = ds„ 
tot s 


entropy, the relation that exists is 
_ i /!« dp tot\ 


ds. 


tot 


dH, 


tot 


L tot 


Jp. 


tot, 


(B2) 


(B3) 


where T^ and are an average between stagnation conditions from one station 

to the next over which Ap^ is being determined. This change in entropy is equal to the 
sum of the external entropy change ds g due to external heat transfer and the internal ir- 
reversible entropy change ds^ due to frictional losses. 

The external entropy change is calculated from 

ds = -^L (B4) 

wT s 


where by definition 


dq _ 

w 

is the increase of coolant total enthalpy between stations. 



(B5) 
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The internal entropy change is determined from the pressure-energy equation of ref- 
erence 11 


P s \ dV b 

d Pg + + Jp s t s d Sint = 0 

o 


and the differential momentum equation 


2p s Vrf M p V, dV b 

dp c + + _L _° — ° = o 

b 


gD 


H 


g 


from the same reference. 

Subtracting equation (B7) from equation (B6) yields 


Jp s T s ds int = 


2p s V^f At 


gD 


H 


which results in 


ds int = 


2V b f At 
J ^ D H T s 


where T g is an average value over the increment AH. 
Using these results in equation (B3) further yields 


L tot 


Kot " 


d P, 


’tot' 


Jp 


dH tot , 2v b f Af 


toty 


JgD H T 


s 


which reduces to 


dp tot “ Jp tot dH tot l 1 ■ 


T toA 2V b f Mp tot T tot 
T s / gD H T s 


(B6) 


(B7) 


(B8) 


(B9) 


(BIO) 


(Bll) 


where all values of p tot , D H , T tot , T g , and V b are averages over the increment AH. 
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Substituting the equation for friction pressure drop 


ip fr = 


2G av* " 


_ ^s, a v 
p s, av D H, av g gD 


2p « ^jay* A£ 


H, av 


(B12) 


into equation (Bll) yields 


Ap tot " Jp tot, av AH tot f 1 ' * 0t ’ aV ) " Ap fr(““ tj 

\ A s, av/ \ p s, 


p tot, av^tot, av^ 


T 

av s, av 


(B13) 


By establishing a value of p^ at the first station, as discussed in the text of this report, 
successive values of p^ are determined by use of equation (B13). 
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APPENDIX C 


COOLANT-PASSAGE STRESS ANALYSIS 
by Rene E. Chambellan 

To perform the desired analysis without creating a statically indeterminate problem, 
the following simplifying assumptions were made: 

(1) All materials remain elastic. (In an actual situation it is anticipated that local- 
ized regions would yield, due particularly to high thermal stresses. ) 

(2) Bending of the pressure shell is not considered. (This is a very localized effect 
in the throat and knuckle regions. ) 

(3) The shell is much stiffer than the tubes and the tubes therefore will deform as the 
shell, thus offering negligible resistance to shell deformations. 

(4) The tube radius is constant along the length of the pressure shell cone considered. 

(5) The material is the same for the shell and tubes. 

(6) The complete nozzle is uniformly cooled down to the bulk coolant temperature at 
startup. 

(7) No residual stresses remain from the fabrication processes. 

(8) Nuclear heating of the pressure shell is neglected. 

In addition, the nozzle pressure shell is divided into conical sections compatible with 
the heat -transfer and fluid-flow sections. The analysis is repeated at each section and 
the results that influence the next section are carried over. 


Shell 


The hot-gas pressure forces produce the membrane forces per unit length NS and 
N, shown in figure 23, in the shell structure. These forces are defined by equations 
(Cl) and (C2) which are revised equations from reference 13 (p. 112): 


NSj = *2 * NS l(f) COB(, ’2 - '"I 1 


(Cl) 


N 2 = Pg 2 s tan ^2 ( C2 ) 

The last term in equation (Cl) is the reaction of the previous section acting on the section 
under analysis. 
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Pressure shell 



Figure 23. - Nozzle pressure shell membrane forces. Cone slant height from previous station, L; cone 
slant height from present location, S; membrane force in longitudinal direction, NS; membrane force 
in circumferential direction, N; shell thickness, w 


The resulting shell strains are 


%, 2 - £-~— < NS 2 - W 

^shsh, 2 


(C3) 


e h, 2 - <N 2 - V'V 

^sl^sh, 2 


(C4) 


No thermal strains appear in the shell with respect to the tube wall because the tubes and 
the shell are cooled down to the bulk coolant temperature uniformly, and then the tubes 
are heated relative to the shell with the nuclear heating of the pressure shell being ne- 
glected. 


Tubes 


The tubes are first considered to act as if they were not restrained, and then the 
shell restraint with the resulting stresses is added. The tubes are forced to match the 
strain that the shell restraint dictates. 

The pressure difference across the tube crown produces a membrane stress of 


a 


P 


(p £ - p g ) R t 


w 


(C5) 


The temperature gradient, assumed linear, through the tube crown produces a biaxial 
stress at the surfaces as follows: 
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O^rp "f" 


¥ (T », g - T .,i | 

2(1 - y 


(C6) 


where compression exists at the outer surface and tension at the inner surface as denoted 
by the minus and plus signs, respectively (see ref. 14, p. 337). 

The tube axial and tangential strains due 
to pressure and heating (see fig. 24) are 



Figure 24. - Stress model for tube axial and tangential strain due to 
pressure and temperature. 


The stress produced by this strain is 


y p(Pj> - Pg )R t 




+ a av^met ” 


(Pf - Pg) R t 


(C7) 

(C8) 


e 0 + “av^ T met T J^ 

E t t w 

The crown tube temperature T met is used 
to obtain the average tube growth. 

By considering the net restraint on the 
tube produced by the shell, the longitudinal 
strain in the tube is 


6 long " e M 6 A 


(C9) 


a long E t e long 

The total longitudinal stress is composed of the combination of this stress plus the bend- 
ing stress due to the temperature gradient through the tube wall 


a long, com a long + a T E t e long 


E t^ T w,g T w,l> 


2d - , p ) 


The tangential deflection of the tube produced by the shell restraint is 

5 0 = [ e 6 " ^p^M " e A^] 2R t ■ ^ D sh " 2R P e h 


(CIO) 
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From reference 14 (p. 156) case 1, the required forces and resulting moment acting 
on the tube to produce this deflection are 


W = 


Ww 


1.788 r; 


K 1 - $ 


(Cll) 


and 


M max = °' 3183 WR t (C12) 

with the maximum moment occurring at the tube crown (see fig. 25). Also, from refer- 
ence 14, the resulting circumferential bending stress at the tube crown is 

“bend " 16 <C13) 

‘w 

which is tension at the outer surface and compression at the inner surface. 

The membrane stress at the tube crown from figure 25 is 



(a) Bending moment. 


Wall thickness, t, 


Deflection force, W 


Gas pressure, p g 



-( p J- p g) R t 

Tube radius, 


Liquid pressure, P| 

I 

HH p J- p g) R t 


CD-8585 


-W 12 

(b) Membrane forces. 

Figure 25. - Deflection force resulting in bending and membrane stresses. 


a M = a p + a d = 


(P £ - P g )R t 


w 


W 

2t w 


or 


a M = 


_ (P£ - Pg) R t 


w 


Ww 


3. 576 R 


?(‘ - 4 ) 


(C14) 


The total tangential stress at the 
tube crown is then the sum of the bending 
stress due to the tube deflection, the 
thermal stress due to the temperature 
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gradient through the tube crown, and the membrane stress due to pressure and tube de- 
flection: 



(C15) 

(C16) 


The material properties E^, E^., 
at the appropriate temperatures. E^ 


v , and a used in these equations are required 
and Vp are supplied as input and should be val- 


ues for the shell and both shell and tube, respectively, at an average temperature in the 


nozzle. Ej. is evaluated at the average metal -wall temperature, and a is evaluated at 
the average temperature of the associated temperature differences. 
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APPENDIX D 


DETAILS OF DESIGN AND EVALUATION PROGRAMS 

The two programs presented in this report are written in FORTRAN IV for use on an 
IBM 7094 Model 2 utilizing the IBSYS (International Business System) system. The pro- 
grams have an average run time of 0. 5 to 2. 0 minutes, largely dependent on the number 
of sections into which the nozzle is divided. The input and output formats utilized by both 
programs are the same except for minor variations as explained later in this section. 

The basic equations utilized in both programs are the same with internal variations 
in procedure and variable names. However, the major variables used throughout the pro- 
grams are consistent between both programs. 


Data Input 

The basic means of data input is through the use of the NAMELIST statement using 
the title STUFF. The variable names and definitions are shown in tables IV and V for the 
array and single valued variables, respectively. One significant factor that must be ob- 
served in the input of data for a variable in an array, is that the first element of the array 
must contain a zero. The only exception to this rule is for the input of the variable X 
into a double array, here no leading zero is used. The first element of the variable array 
is used as a computational location and the other elements are moved into this location as 
they are needed for computations. 

The material properties k, a, E, and a y as a function of temperature are input in 
the form of coefficients to a polynomial ranging to the ninth degree and located in the 
KMET subroutine. The coefficients of the polynomials were obtained by a least-squares 
method of curve fitting applied to data from reference 15 and to unpublished data. The 
resulting coefficients are those shown in table VI with the resulting curves being those 
previously shown in figure 6 (p. 20). For the coefficients presented in table VI, the tem- 
perature range of applicability is from 300° to 2200° R. These limits are input in the 
form of TEMPI and TEMP2, respectively, with TEMP3 being an upper limit on the tem- 
perature for the thermal conductivity of a coating material. In the event one of these 
limits is exceeded, the limiting temperature value is used to determine the properties. 

The hydrogen fluid thermodynamic and transport properties used in both programs 
are obtained from a FORTRAN IV version of the subroutine described in reference 16. 

The properties subroutine is called from a library within the Lewis Computer System and 
would have to be added as a subroutine on other systems. The hydrogen composition used 
in the programs is specified as input data in the form of the variable COMP. 
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TABLE IV. - ARRAY VARIABLES FOR NAMELIST STUFF 


FORTRAN 

symbol 


ANWEB 

AREACR 2 

AREAL b 

AX 

COEFG 

COEFL 

CT 

DHG 


PGS 
| PHI 
QGEN 
QRAD 
RC 

ROUGH a 

SHELTH 

TESTC C 


TGS 

TKCER C 

TW 

TWG b 

IX 


Definition 


Number of coolant tubes at station 
Area correction, e Q 
Coolant flow area per tube 
Axial length from datum plane to station 
Hot-gas correlation coefficient 
Coolant correlation coefficient 
Convergence tolerance used 
Engine diameter to gas side of metal wall at 
tube crown 
Hot-gas static pressure at station 
Angle gas wall makes with nozzle centerline 
Nuclear heating rates in shell 
Radiation heat flux entering gas wall 
Radius of gas-wall curvature in axial direction 
Tube-wall roughness on coolant side 
Average pressure shell wall thickness 
Coating test, negative for no coating, positive 
for given coating, and zero for thermal 
barrier 

Hot-gas static temperature at station 
Coating thickness or thermal resistance 
Metal-wall thickness 
Hot-gas-side metal-wall temperature 
Coefficients for material properties polynomial 


Units 


sq in. 
in. 


see table VII 
in. 

psia 

deg 

Btu/(in. 3 )(sec) 

Btu/(in. 2 )(sec) 

in. 

in. 

in. 


in. or (in. 2 )(sec)(°R)/Btu 


in. 

R 


a AREACR and ROUGH are used as a single value in the design program, and a value 
is specified at each station in the evaluation program. 
b AREAL is used only in the evaluation program input, and TWG is used only in the 
design program input. 

c TESTC and TKCER are used only as input to the evaluation program. 
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TABLE V. - SINGLE VALUED VARIABLES FOR NAMELIST STUFF 


FORTRAN 

Definition 

Units 

symbol 



AKR a 

Estimate used in determining recovery factor 


AMU 

Average Poisson ratio of material 


BRAZ 

Thickness of braze material 

in. 

Cl b 

Option on use of Cl correlation 


C2 b 

Option on use of C2 correlation 


C3 b 

Option on use of C3 correlation 


COMP 

Fraction of ortho -hydro gen present in composition 


DEN 

Average density of tube and shell material 

lb/in. 3 

DIAM 

Hot-gas-side nozzle-throat diameter 

in. 

ES 

Average Youngs modulus of shell at its operating temperature 

psi 

Fl c 

Estimate used in determining friction factor 


NUMB 

Number of sections nozzle is divided into 


PLIN 

Coolant inlet static pressure 

psia 

QEST d ’ e 

Estimate used in determining convection heat flux 

Btu/(in. ^)(sec) 

TCER d,f 

Initial estimate of average coating temperature 

°R 

TEMPI 

Lower limit of temperature range on material properties 

°R 

TEMP3 d 

Upper limit of temperature range on material properties 

°R 

TEMP3 

Upper limit of temperature range on coating properties 

°R 

TLIN 

Coolant inlet total temperature 

°R 

TMETl d,g 

Initial estimate of average metal -wall temperatures 

°R 

TWGE d ’ b 

Initial estimate of hot-gas-side metal-wall temperature 

°R 

WG 

Hot -gas mass flow rate 

lb/sec 

WL 

Coolant mass flow rate 

lb/sec 

YM b 

Mach number limit 



a Suggested initial estimate for use in evaluation program, 0.9. 
b Used only in design program. 

c Suggested initial estimate for use in evaluation program, 0. 003. 
d Used only in evaluation program. 
e Suggested initial estimate, 1.0. 

Suggested initial estimate, 1500.0. 

^Suggested initial estimate, 750.0. 

^Suggested initial estimate, 1000.0. 


TABLE VI. - POLYNOMIAL COEFFICIENTS FOR 


PROPERTY CURVES OF 347 STAINLESS STEEL 


KMET a 

Thermal 

conductivity 

Yield 

strength 

Young’s 

modulus 

Thermal 

expansion 

A 

1.456672 

37.72398 

20.60063 

5.05081 

B 

1. 131624 

-9.03007 

62.24962 

25.84439 

C 

-1.341449 

-5.64351 

-168.24084 

-71.93478 

D 

1.912256 

3.49807 

214.93276 

111.25141 

E 

-.951701 

-.90454 

-155.08070 

-98.19083 

F 

. 155014 

0 

62.56724 

49.60798 

G 

0 

0 

-12.96823 

-13.37053 

H 

0 

0 

1.03951 

1.49251 

I 

0 

0 

0 

0 

J 

0 

0 

0 

0 

Scale fac- 
tor, SF 

l.OxlO" 4 

1.0X10 3 

l.OxlO 6 

l.OxlO" 6 


a KMET = (A + BT + CT 2 + . . . + JT 9 )SP) where T is 
the temperature in °R divided by 1000. This division 
and the scaling factor are taken care of internally 
within the programs. 

The major input data difference between the two programs is that the design program 
uses the hot-gas-side metal-wall temperature as input, and the evaluation program uses 
the exact coolant passage area as input. Other differences are that the design program 
accepts only single values of area correction e c and tube roughness as input while the 
evaluation program accepts values at each individual station; the means of obtaining dif- 
ferent options as to the use of the Cp C 2 , Cg, and rough or smooth tube conditions 
varies as shown in table VII; the evaluation has the flexibility of handling both the radial 
and flat plate cases of heat conduction across the tube crown and a choice must be speci- 
fied as input data. One other difference exists in the input of data, that difference being 
that the evaluation program utilizes realistic initial estimates of certain parameters at 
the first station, as shown at the end of table V. 


Data Output 

The basic output format is shown at the end of this appendix for case 1 (p. 65). The 
input data is printed out for the purpose of maintaining run records and checking. The 
definitions of the output variables that are not self-explanatory are shown in table VUE. 

In addition to the output shown additional write-outs are possible to note the use of a tube 
splice, in the event the materials properties temperature range is exceeded, in the event 
the coolant flow chokes or an iteration loop does not converge. 
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TABLE Vn. - PROGRAM OPTIONS 


LESTC1 

LESTC2 

LESTC3 

LESTC4 

LESTC5 

Cl 

C2 

C3 

ROUGH 


Evaluation logicals 

True for use of 

True for use of C 3 

True for use of C 3 

True for assumption of smooth tubes 

True for assumption of radial heat conduction 

Design program reads under NAMELIST 

0.0 or 1.0; 1. 0 for use of Cj 
0. 0 or 1. 0; 1. 0 for use of Cg 
0. 0 or 1. 0; 1. 0 for use of C 3 
0. 0 or value; 0. 0 for use of smooth tubes 


TABLE VIE. - OUTPUT PARAMETERS 


Output headings 

Definition 

Q/A OUT 

2 

Heat flux entering outside surface of tube crown, Btu/(in. )(sec) 

Q/A IN 

2 

Heat flux leaving inside surface of tube crown, Btu/(in. )(sec) 

T. REF. 

Temperature utilized in determining hot-gas properties, °R 

T.COAT.G 

Temperature on hot-gas side of coating or thermal barrier, °R 

T. MET. G 

Temperature on hot- gas side of metal wall, °R 

T. MET. L 

Temperature on coolant side of metal wall, °R 

Q/SECTION 

Total heat added to coolant between stations, Btu/sec 

Q/A NUC 

Heat flux passing through shell wall, Btu/(in. )(sec) 

HYD.DIA 

Coolant passage hydraulic diameter, in. 

WALL AREA 

2 

Total heat -transfer surface area, in. 

YIELD STRENGTH 

Yield strength at average metal-wall temperature, psi 

COE F. OF LIN. EXP 

Coefficient of linear expansion at average temperature between 


metal wall and shell, in. /in. 

YOUNG’S MODULUS 

Young’s modulus at average metal-wall temperature, psi 






















Figure 26. - Design program flow diagram. INumbers correspond to FORTRAN statement numbers. ) 
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Figure 27 . - Evaluation program flow diagram. {Numbers correspond to FORTRAN statement numbers. ) 



















Program Details 


The main programs utilize the basic equations set forth in the text of this report. 

The FORTRAN IV coding of the design (pp. 70 to 79) and evaluation (pp. 80 to 89) pro- 
grams is included. The basic flow diagrams of both the design and evaluation program 
are shown in figures 26 and 27, respectively. The iterative procedures contained in these 
programs require specified convergence tolerances for solution. These convergence tol- 
erances are defined in table EX and are used as fixed internal values as shown unless 
changed by specified input values. Internal limits are also specified as to the maximum 
number of iterations required for convergence. The design program accepts the final 
value if convergence has not been reached and proceeds with the calculations. The eval- 
uation program, however, stops and prints an appropriate write-out and the answers cal- 
culated to that point. The evaluation program also employs a straight line iteration sub- 
routine called ITRAT to determine the next estimate during some of the iteration. 


TABLE IX. - INTERNAL CONVERGENCE TOLERANCES UNLESS CHANGED 


FORTRAN 

symbol 

Definition 

Value 

Units 

CT( 1) 

Used to solve equations (19) and (20) for friction factors 

3.0xl(T 6 


CT(2) 

Used to determine the coolant -passage height and also 
width in design program 

2.0X10' 4 

in. 

CT(3) 

Used to determine average metal temperature, average 
coating temperature in evaluation program, and hot- 
gas film temperature in design program 

1 

°R 

CT(4) 

Used to determine static coolant pressure and also total 
coolant pressure in evaluation program 

0. 1 

psi 

CT( 5) 

Used to determine hot-gas recovery factor 

l.OxlO " 3 


CT(6) 

Used to determine static coolant enthalpy 

l.OxlO " 4 

Percent 

difference 

CT(7) 

Used to determine heat-flux balance across tube wall 

l.OxlO " 3 

l.OxlO " 3 

a Btu/(in.^)(sec) 

^Percent 

difference 

CT(8) 

Used to determine total coolant temperature in design 
program and gas -side metal -wall temperature in 
evaluation program 

0. 1 

°R 

CT(9) 

Used to determine coolant state at first section in 
evaluation program 

l.OxlO -4 

Percent 

difference 


a E valuation program. 
b Design program. 
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The programs are set up to handle counterflow of the coolant and hot gas, but could 
be revised or used for parallel flow of the coolant and hot gas. The revision would re- 
quire a minor variation in the output of data and a change in the sequence of stress calcu- 
lations to allow the integration of the hot-gas pressure forces acting on the shell to pro- 
ceed from the hot-gas exit end back toward the chamber. Because the heat -transfer and 
fluid-flow calculations are not interrelated to the stress calculations, the program can be 
used for parallel flow if the resulting stress values are not considered. 

As mentioned previously in this report, both programs have the provision for incor- 
porating a tube splice in the form of a bifuraction. For the purpose of analysis, a station 
is chosen at the point of the tube splice. The number of tubes specified at the transition 
station is that number of tubes which exists between the preceding station and the transi- 
tion station (see fig. 28). In addition, the length between stations, where the transition 
occurs, should be as small as possible while including all the region of uncertain tube 
geometry. 

Two internal program variations that exist are associated with the average density 
used to calculate the coolant inlet total pressure and the method of handling the calcula- 


Tube ends formed to produce 
no variation in tube wall 



Section D-D 


tion of pressure drop over the transition 
increment for the situation where a tube 
splice exists. The design program uses 
the static density instead of an average 
density for calculating total pressure at 
the first station (eq. (24)). The pressure- 
drop calculations use the areas at both the 
preceding and present stations. Over the 


1 

r*~ A 

— B / 
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Station i + 1 


1 

i 

> 

- \ 

Station i 

i 

i 

\ 

\ 

1 

Shell side 

^A 



L c ' 



Section E-E 



This flow area is 
modified for pressu re- 



section A-A Section B-B Section C-C 

CD-8582 

Figure 28, - Tube splice transition section showing geometry before, at, 
and after bifurcation. 


spliced-tube transition section, the flow 
from two tubes enters one tube or vice 
versa depending on where the splice is 
located. Therefore, the change in flow 
area from one station to the next is not 
actually the change in area from that of 
one tube before the splice to that of one 
tube after the splice. A revision of the 
previous area is made in each program to 
account for this change in the number of 
tubes. The design program assumes that, 
at the splice joint, the area and radius of 
one small U-tube equals one -half the area 
and radius of one large U-tube, and dimen- 
sional revisions for the pressure -drop 
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calculations are made accordingly. The evaluation program assumes that the area, ra- 
dius, and geometrical configuration are as shown in figure 28. This is more nearly the 
actual condition where two small tube ends would be formed into a back-to-back symmet- 
rical configuration to fit into one large U-tube. 

The evaluation program has one option for restart, that option being on the conditions 
which could result in choked coolant flow. This option is exercised when the static cool- 
ant pressure drops to 0. 3 of the inlet coolant pressure. The option causes a written mes- 
sage of output to be printed, increases the coolant inlet pressure by 50 pounds per square 
inch, and restarts the evaluation run at the first station. The design program has no pro- 
gram restarts, but does contain other coolant and geometry limits as explained in the 
section Basic Analytical Procedure. 

As explained in the section Design Procedure, care should be exercised in subdividing 
the nozzle into sections, because excessively large sections can lead to errors in the pres- 
sure calculations, trouble in excessive iterating of the pressure calculation, and the over- 
looking of a nozzle hot spot. Also, using large section lengths in the convergent section, 
where the slope of the heat flux curve is steep, can lead to a problem in the evaluation 
program. The problem occurs in the heat -flux iteration that uses the value of heat flux 
from the last station as an estimate to start calculations at the next station. A sharp drop 
in heat flux results in the use of an excessively large heat -flux estimate that produces un- 
realistic temperatures in the program. Traps are present in the program to work this 
condition out, but if conditions become excessive the iteration limits could be exceeded. 

The programs are set up so that cases can be executed in sequencial order by simply 
stacking the input data back to back. Variations of input data for the same nozzle can be 
obtained by simply reading in only the input values that vary and the logicals in the evalu- 
ation program. Each program uses the last value of friction factor and recovery factor 
calculated for a given case for input to the succeeding case unless their input values are 
also repeated. These calculated values are just as adequate as the input estimates used 
initially and need not be repeated. For the design program, if the input values for gas- - 
wall temperature are not achieved and new values are calculated, these values will be 
used as input to the succeeding case unless the original values are again specified as in- 
put. 

Cases 3, 6, and 7 were run by using the reference 4 correlation on the hot-gas side 
with internal program revisions in statements 110, 112, and 122. The existing state- 
ments in the evaluation program were replaced by 

110 TGF = (TGE + TGCER)/2. 0 
112 TGF = (TGE + TWGE)/2.0 
122 PRGEX = PRGF **0. 4 
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BASIC OUTPUT FORMAT 


FLUID PROPERTY LIBRARY WITH A COMPOSITION OF 96. B PERCENT PARA-, AND 3.2 PERCENT ORTHO— HYOROGEN, Ll/I/61 


JCONST 

CT 

3.0000000E-06, 

2. OOOOOOOE— 04- , 

I.OOOOOOOE 00, 

1. OOOOOOOE-Ol, 

I.OOOOOOOE- 03, 

1 . OOOOOOOE— 04 , 

X 

1. 00000 00E-03, 
1.4566720E 00, 

l. 0000000E-01 , 
1.1316240E 00, 

1.0O00000E-04, 
-1.3414490E 00, 

1.9122560E 00, 

-9, 51701 G0E-01 , 

l. 5501400E— 01 , 


0. , 

0* , 

0 . , 

0. , 

3.7724000E 01, 

-9.0301000E 00, 


-5.6435000E 00, 

3. 4981 00 OE 00, 

-9.0450000E-01, 

0 . , 

0 « , 

0* , 


0* i 

0. , 

2.0500600E 01, 

6.2249600E 01, 

— 1 . 6824080E 02, 

2. 1493280E 02, 


- 1 .55 08070E 02, 

6. 2567200E 01, 

— 1 . 2968200E 01, 

1.0395000E 00, 

0. , 

0. , 


5.05 08000E 00, 

2* 5844400E 01, 

-7.19348006 01, 

1.1125140E 02, 

—9. 8190799E 01, 

4. 9607999E 01, 


-1.3370500E 01, 

1. 492 5000E 00, 

0. , 

0. , 

0. , 

0. , 


0. , 

0. , 

0. 

0. , 

o. , 

0. 

TESTC * 

0. , 

0 • , 

0. 

-5. OOOOOOOE-Ol , 

-5.00000006-01, 

-5.0000000 E-01, 

-5. OOOOOOOE-Ol, 

-5. OOOOOOOE-Ol* 


-5. OOOOOOOE-Ol, 

-5. OOOOOOOE— 01 , 

-5.0000003 E-Ol, 

-5 .OOOOOOOE-Ol, 

-5. OOOOOOOE-Ol, 

-5. OOOOOOOE-Ol , 


-5. OOOOOOOE-Ol, 

-5. OOOOOOOE-Ol , 

-5. OOOOOOOE-Ol, 

-5 .OOOOOOOE-Ol, 

-5. OOOOOOOE-Ol, 

-5. OOOOOOOE-Ol, 


-5 .OOOOOOOE- 01 , 

-5. OOOOOOOE-Ol , 

-5. OOOOOOOE-Ol , 

-5. OOOOOOOE-Ol, 

-5. OOOOOOOE-Ol, 

-5* OOOOOOOE-Ol , 


-5.00000 00E-01, 

-5. OOOOOOOE -01, 

0. ' 

0 . , 

0. , 

0. , 


0 • , 

0. 

0. 

0. , 

0. , 

0. , 


0. , 

0* , 

0. , 

0. , 

0. , 

0. 


0. , 

0. , 

0. 

0 . , 

0. , 

0. , 


0. , 

0. 






TEMPI = 3.0000000E 02, TEMP 2 » 2.2000000E 03, TEMP3 * 0. 


QEST = i.OOOOOOOE 00* TCER * 1.5000000E 03, 


TWGE = I.OOOOOOOE 03, TMETl = 7.5000000E 02, FI 


3.0000000E-03, AKR 


8.9999999E-01, 


* END 


SEC.NO. 

AXIAL LENGTH 

ENG .01 AM. 

NO. TUBES 

ANGLE PHI 

TUBE THICK. 

AREA PER .TUBE 

STAT.GAS TEMP 

STAT.GAS PRESS 


1 IN.I 

CIN. 1 


(DEG.) 

( IN.) 

( SQ.IN.) 

(DEG. R 1 

fPSI A) 

1 

0 . 

12.1260 

90.0 

15.000 

0.0120 

0.0980 

1225.0 

5.90 

2 

2.0000 

11.0540 

90.0 

15.000 

0.0120 

0.0815 

1325.0 

7.80 

3 

4.0000 

9.9820 

90.0 

15.000 

0.0120 

0.0671 

1445.0 

10.60 

4 

6.0000 

8.9100 

90.0 

15.000 

0.0120 

0.0539 

1600.0 

15.90 

5 

8.0000 

7. 8380 

90.0 

15.000 

0.0120 

0.0419 

1795.0 

24.50 

6 

10.0000 

6.7680 

90.0 

15.000 

0.0120 

0.0317 

2040.0 

40.00 

7 

11.0000 

6.2300 

90.0 

15.000 

0.0120 

0.0268 

2200.0 

50.00 

8 

12.0000 

5.6940 

90.0 

15.000 

0.0120 

0.0227 

2400.0 

68.00 

9 

13.0000 

5.1580 

90.0 

15.000 

0.0120 

0.0187 

2635.0 

96.00 

10 

13.5000 

4. 8900 

90.0 

15.000 

0.0120 

0.0167 

2800.0 

120.00 

11 

13.9000 

4.6760 

90.0 

15.000 

0.0120 

0.0153 

2940.0 

148.00 

12 

14.2240 

4.5040 

90.0 

15.000 

0.0120 

0.0142 

3100.0 

184.00 

13 

14.7000 

4.3300 

90.0 

5.750 

0.0120 

0.0130 

3325.0 

240.00 

14 

15.0000 

4.3000 

90.0 

0. 

0.0120 

0.0125 

3458.0 

287.00 

15 

15.3000 

4. 3300 

90.0 

5.750 

0.0120 

0.0125 

3598.0 

323.00 

16 

15.8000 

4. 5180 

90.0 

15.480 

0.0120 

0.0128 

3643.0 

353.00 

17 

16.4000 

4.9920 

90.0 

27.800 

0.0120 

0.0142 

3852.0 

448.00 

18 

17.1210 

6.0580 

90.0 

45.000 

0.0120 

0.0182 

3942.0 

500.00 

19 

18.4000 

8. 6160 

90.0 

45.000 

0.0120 

0.0335 

3987.0 

522.00 

20 

20.0760 

11.9680 

90.0 

45 .030 

0.0120 

0.0620 

3996.0 

528.00 

21 

21.8000 

14.2160 

90.0 

21.300 

0.0120 

0.0900 

3998.0 

528.00 

22 

23.6120 

14. 8960 

90.0 

0 . 

0.0120 

0.1000 

3998.0 

530.00 

23 

25.6000 

14.6960 

90.0 

0 . 

0.0120 

0.1000 

3999.0 

530.00 

24 

28.0000 

14.8960 

90.0 

0 . 

0.0120 

0.1000 

4000. 0 

530.00 

25 

31.0000 

14. 8960 

90.0 

0 . 

0.0120 

0.1000 

4000. 0 

530.00 
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SEC. NO. 

GAS COEF. 

COOL. COEF. 

AREA CORR. 

ROUGHNESS 

(IN.) 

RAO. OF CJRV 
( IN.) 

Q RAD. 

0TU.PER. 
SQ . IN— SEC 

NUC.HEAT 

8TU.PER. 
CU. IN- SEC 

SHELL. THICK. 
(IN.) 

COAT. THK. OR RES. 

UN) OR SO. IN-SEC-R 
PER.BTU. 

1 

0.0260 

0.0210 

0.8000 

0.600 E-04 

0.10006 

08 

0 . 

0 . 

0.750 

0. 

2 

0.0247 

0.0210 

0. 8000 

0.600E-04 

0.1000E 

08 

0. 

0. 

0.750 

0. 

3 

0.0240 

0.0210 

0 . 8000 

0.600E-04 

O.IOOOE 

08 

0. 

0. 

0.750 

0. 

4 

0.0230 

0.0210 

0.8000 

0.600E-04 

0.1000E 

08 

0. 

0 . 

0.750 

0. 

5 

0.0220 

0.0210 

0. 8000 

0.600E-04 

0.1000 E 

08 

0 . 

0. 

0.750 

0. 

6 

0.0209 

0.0210 

0. 8000 

0.600E-04 

0.10006 

08 

0 . 

0 . 

0.750 

0 . 

7 

0.0203 

0.0210 

0. 8000 

0.600E-04 

0. 1000 E 

08 

0. 

0. 

0.750 

0 . 

8 

0. 01 98 

0.0210 

0. 8000 

0.6006-04 

0.1000E 

08 

0 . 

0 . 

0.750 

0. 

9 

0.0193 

0.0210 

0. 8000 

0.600E-04 

o.ioOoe 

08 

0 . 

0 . 

0.750 

0. 

10 

0.0191 

0.0210 

0.8000 

0.6 00 E-04 

0 .1000 E 

08 

0. 

0 . 

0.750 

0 . 

11 

0.0192 

0.0210 

0.8000 

0 .6006—04 

O.IOOOE 

08 

0. 

0 . 

0.750 

0. 

12 

0.0194 

0.0210 

0. 8000 

0.600 E— 04 

O.IOOOE 

08 

0. 

0 . 

0.750 

0. 

13 

0.0200 

0.0210 

0. 8000 

0.600E-04 

0.3000E 

01 

0 . 

0. 

0.750 

0. 

14 

0.0202 

0.0210 

0.8000 

0.600E— 04 

0.3000E 

01 

0. 

0 . 

0.750 

0 . 

15 

0.0204 

0.0210 

0.8000 

0. 500 E-04 

0 .3000 6 

01 

0. 

0 . 

0.750 

0. 

16 

0.0207 

0.0210 

0. 8000 

0.600E-04 

0 • 3000 £ 

01 

0. 

0 . 

0.750 

0. 

17 

0.0242 

0.0210 

0.8000 

0 .600 6—04 

0.3000 E 

01 

0. 

0 . 

0.750 

0 . 

18 

0.0288 

0.0210 

0. 8000 

0 ,600 E-04 

0 • 3000 E 

01 

0. 

0 . 

0.750 

0 . 

19 

0.0316 

0.0210 

0. 8000 

0.600F-04 

0,10006 

08 

0. 

0 . 

0.750 

0. 

20 

0.0321 

0.0210 

0. 8000 

0 .6006-04 

O.IOOOE 

08 

0. 

0 . 

0.750 

0. 

21 

0.0322 

0.0210 

0. 8000 

0.5006-04 

-0.5000 E 

01 

0. 

0 . 

0.750 

0 . 

22 

0.0323 

0.0210 

0. 8000 

0.6006-04 

-0.50006 

01 

0. 

0. 

0.750 

0 . 

23 

0.0323 

0.0210 

0. 8000 

0.6006-04 

-0.10006 

08 

0. 

0 . 

0.750 

0. 

24 

0.0323 

0.0210 

0.8000 

0.500E-04 

-O.IOOOE 

08 

0. 

0. 

0.750 

0 . 

25 

0.0323 

0.0210 

0. 8000 

0.6006-04 

-O.IOOOE 

08 

0. 

0. 

0.750 

0 . 


AVG.HET AL DENS IT Y= 
PROPELLANT 

0.280 LB S/C U. I N SHELL YOUNG* S MODULUS* 0.300E 

FLOW RATE 16.6650 LBS/SEC 

08 PSI POISSON'S RATIO* 0.290 

COOLANT FLOW RATE 16.6650 

BRAZE THICK- 0.0100 (IN. 
LBS/SEC 

CHAMBER 

TEMPERATURE 

4000.0 DEG. R 

ENTRANCE COOLANT TEMP. 56.00 

OEG.R 

CHAMBER 

PRESSURE 

530.0 PSIA 

ENTRANCE COOLANT PRESS 1150.0 

PSIA 




THROAT DIAMETER 
TOTAL NOZZLE LENGTH 
CONVERGENT AREA RATIO A/A* 
DIVERGENT AREA RATIO A/ A* 

4.3000 I IN • 1 
31.0000 ( IN.) 

12.0006 
7.9524 





C-l rS NOT USED 

C-2 IS NOT USEO 

C-3 IS NOT USED 

TUBE IS ASSUMED SMOOTH 

ASSUMED RADIAL HEAT CONDUCTION 



S EC .NO . 

Q/A OUT 
( 8 TU. /SO. IN 

Q/A IN 
.-SEC. ) 

EFF.GAS T 
IDEG.R) 

T .REF. 
IDEG.R) 

T.COAT.G 

(DEG.R) 

t.met.g 

(DEG.R) 

T .MET .L 
( DEG.R) 

LIQ.T. 
TOTAL 
( DEG.R ) 

LIQ.P. 

STATIC 

(PSIA) 

LIQ.VEL 

FT/SEC 

MACH NO 

1 

1.62 

1.72 

3644. 

1880. 

0 . 

1471. 

1401. 

56.0 

1150. 

63. 

0.016 

2 

1.93 

2.05 

3658. 

1912. 

0 . 

1472. 

1389. 

59.6 

1149. 

78. 

0.020 

3 

2.34 

2.52 

3676. 

1970. 

0 . 

1514. 

1414. 

63.3 

1147. 

98. 

0.026 

4 

2.91 

3.15 

3699. 

2038. 

0 . 

1553. 

1429. 

67.2 

1144. 

126. 

0.035 

5 

3.72 

4.08 

3726. 

2123. 

0 . 

1601. 

1444. 

71.3 

1138. 

170. 

0.049 

6 

4.87 

5.43 

3759. 

2233. 

3. 

1669. 

1465. 

75.7 

1126. 

237. 

0.072 

7 

5.69 

6.40 

3780. 

2301. 

0 . 

1707. 

1469. 

78.0 

1114. 

289. 

0.091 

8 

6.68 

7.61 

3 806. 

2400. 

0 . 

1780. 

1504. 

80.4 

1097. 

353. 

0.114 

9 

7.96 

9.22 

3836. 

2512. 

0 . 

1861. 

1535. 

82.9 

1068. 

446. 

0.149 

10 

8.79 

10.28 

3856. 

2591. 

0 . 

1917. 

1559. 

84.1 

1044. 

511. 

0.174 

11 

9.54 

11.24 

3873. 

2669. 

0 . 

1987. 

1604. 

85.2 

1020. 

569. 

0.198 

12 

10.24 

12.15 

3 893. 

2750. 

0 . 

2067. 

1660. 

86.1 

996. 

625. 

0.221 

13 

11.17 

13.36 

3920. 

2883. 

0 . 

2180. 

1744. 

87.4 

958. 

703. 

0.256 

14 

11. 56 

13.85 

393 6. 

2940. 

0 . 

2211. 

1762. 

88.2 

935. 

745. 

0.276 

15 

11.71 

14.01 

3953. 

2 996. 

0 . 

2239. 

1786. 

89.0 

925. 

758. 

0.284 

16 

11. 55 

13.71 

3958. 

2981. 

0. 

2180. 

17 30. 

90.5 

917. 

762. 

0.290 

17 

11.38 

13.26 

3983. 

3075. 

0 . 

2241. 

1808. 

92.4 

932. 

710. 

0.273 

18 

10.06 

11.38 

3993. 

3 079. 

0 . 

2193. 

1815. 

95.3 

971. 

578. 

0.223 

19 

6.64 

7.22 

3998. 

2996. 

0 . 

2000. 

1748. 

101.6 

1018. 

345. 

0.133 

20 

4. 19 

4.44 

4000. 

2899. 

0 . 

1800. 

1637. 

109.4 

1033. 

210. 

0.081 

21 

3. 19 

3.35 

4000. 

2867. 

0 . 

1736. 

1611. 

115.3 

1037. 

156. 

0.061 

22 

3.01 

3.15 

4000. 

2845. 

0 . 

1692. 

1572. 

120.3 

1037. 

150. 

0.058 

23 

3.09 

3.24 

4000. 

2819. 

0 . 

1638. 

151’. 

126.1 

1036. 

161. 

0.062 

24 

3. 19 

3.34 

4000. 

2789. 

0 . 

1578. 

1445. 

133.4 

1035. 

174. 

0.066 

25 

3.30 

3.46 

4000. 

2 754. 

0 . 

1509. 

1368. 

143.4 

1033. 

192. 

0.072 
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SEC. NO. 

FR I • P .DROP 

MOM. P. DROP 

FRI.FACT 

HGF 

HLF 

Cl 

C 2 

C 3 

O/SECTION 

Q/A NUC 


IPS!) 

(PSD 


(BTU/SQ.1 

N. SEC. DEG. R) 




BTU/SEC 

BTU/SQ. tN.SEC 

1 

0. 

0. 

0. 

0.00075 

0.00128 

1.000 

1.000 

1.000 

0. 

0. 

2 

0.18 

0.96 

0.00296 

0.00088 

0.00155 

1.000 

1.000 

1.000 

169.6 

0. 

3 

0.28 

1.54 

0.00286 

0.00108 

0.00186 

1.000 

1.000 

1.000 

184.9 

0. 

4 

0.47 

2.72 

0.00277 

0.00136 

0.00231 

1.000 

1.000 

1.000 

203.6 

0. 

5 

0.84 

5.23 

0.00267 

0.00175 

0.00297 

1.000 

1.000 

1.000 

226.7 

0. 

6 

1.64 

10.45 

0.0025 6 

0.00233 

0.00390 

1.000 

1.000 

1.000 

254.9 

0. 

7 

1.47 

10.22 

0. 00248 

0.00274 

0.00460 

1.000 

1.000 

1.000 

138.8 

0. 

0 

2.25 

14.89 

0.00242 

0.00330 

0.00535 

1.000 

1.000 

1.000 

148.4 

0. 

9 

3.55 

25.81 

0.0023 6 

0.00403 

0.00635 

1-000 

1.000 

1.000 

159.1 

0. 

10 

2.65 

21.16 

0.00232 

0.00454 

0.00697 

1.000 

1.000 

1.000 

83.9 

0. 

11 

2.75 

20.98 

0.00229 

0.00506 

0.00740 

1.000 

1.000 

1.000 

69.7 

0. 

12 

2.75 

21.69 

0.00226 

0.00561 

0.00772 

1.000 

1.000 

1.000 

58.3 

0. 

13 

4.93 

33.00 

0.00223 

0.00641 

0.00806 

1.000 

1.000 

1.000 

87.4 

0. 

1 4 

3.65 

19. 13 

0.00221 

0.00670 

0.00827 

1.000 

1.000 

1.000 

56.2 

0. 

15 

3.88 

6.05 

0.00220 

0.00683 

0. 00826 

1.000 

1.000 

1.000 

57.6 

0. 

16 

6.52 

1.81 

0.00219 

0.00650 

0.00836 

1.000 

1.000 

1.000 

100. 1 

0. 

17 

7.28 

-22.30 

0. 00219 

0.00654 

0.00773 

1.000 

1.000 

1.000 

134.9 

0. 

18 

6.81 

-46.68 

0.00222 

0.00559 

0.00662 

1.000 

1.000 

1.000 

204.4 

0. 

19 

4.79 

-51.85 

0.00228 

0.00332 

0.00439 

1.000 

1.000 

1.000 

431.0 

0. 

20 

1.61 

-16.38 

0.0023 8 

0.00190 

0.00291 

1.000 

1.000 

1.000 

519.7 

0. 

21 

0.51 

-4.03 

0.00247 

0.00141 

0.00224 

1.000 

1.000 

1.000 

395.0 

0. 

22 

0.29 

-0.37 

0.00251 

0.00130 

0.00217 

1.000 

1.000 

1.000 

333.1 

0. 

23 

0.29 

0.61 

0.00252 

0.00131 

0.00234 

1.000 

1.000 

1.000 

362.6 

0. 

2* 

0.38 

0.78 

0.00252 

0.30132 

0.00255 

1.000 

1.000 

1.000 

450.9 

0. 

25 

0. 52 

1.02 

0.00252 

0.00133 

0.00283 

1.000 

1.000 

1.000 

582.4 

0. 


SEC .NO. 

SHELL RAO 

TUBE RADIUS 

TUBE HEIGHT 

TUBE WIDTH 

HYD.OIA 

WALL AREA 

TOTAL TUBE HEIGHT 


( IN.) 

UN. ) 

(IN.) 

( IN. ) 

(IN.) 

(SO. IN) 

(IN.) 

1 

6.3503 

0.2018 

D • 08 36 

0.4093 

0.3238 

0. 

0.2855 

2 

5.7915 

0.1825 

0.0794 

0.3703 

0.2957 

95.6227 

0.2619 

3 

5.2345 

0.1631 

0.0769 

0.3314 

0.2690 

86.5540 

0.2401 

4 

4.6768 

0. 143 8 

0.0739 

0.2925 

0.2417 

77.4854 

0.2177 

5 

4.1184 

0.1244 

0.0700 

0.2535 

0.2137 

68.4160 

0.1944 

6 

3.5624 

0.1051 

0.0675 

0.2147 

0.1865 

59.3492 

0.1727 

7 

3. 2814 

0.0954 

0.0649 

0.1951 

0.1717 

26.2841 

0.1602 

8 

3.0035 

0.0857 

0.0644 

0.1757 

0.1583 

24.0062 

0.1501 

9 

2.7243 

0.0760 

0.0625 

0.1562 

0.1439 

21.7390 

0.1385 

10 

2.5838 

0.0712 

0. 0606 

0.1464 

0.1360 

10.0193 

0.1318 

11 

2.4725 

0.0673 

0.0600 

0.1386 

0.1302 

7. 6068 

0.1272 

12 

2.3828 

0.0642 

0.0593 

0.1324 

0.1255 

5.8940 

0.1235 

13 

2.2949 

0.0611 

0.0574 

0.1262 

0.1201 

8.1701 

0.1185 

14 

2.2773 

0.0606 

0.0547 

0.1250 

0. 1 1 77 

4.9495 

0.1153 

15 

2.2909 

0.0611 

0.0535 

0.1259 

0.1177 

4.94 95 

0.1146 

16 

2.3790 

0.0644 

0.0481 

0.1321 

0.1189 

8.6043 

0.1126 

17 

2.6065 

0.0727 

0.0401 

0.1480 

0.1244 

11.7589 

0.1129 

18 

3.1218 

0.0913 

0.0279 

0.1839 

0.1383 

19.0667 

0.1192 

19 

4.4237 

0.1371 

0.0145 

0.2748 

0.1824 

51.6203 

0.1516 

20 

6.1336 

0. 1970 

0.0026 

0.3942 

0.2435 

95.9936 

0.1996 

21 

7. 3422 

0.2393 

0.0001 

0 .4786 

0.2925 

107.1460 

0.2394 

22 

7.7124 

0.2522 

0.0002 

0.5044 

0.3084 

107.4853 

0.2524 

23 

7.7124 

0.2522 

0.0002 

0.5044 

0.3004 

118.8072 

0.2524 

24 

7.7124 

0.2522 

0.0002 

0.5044 

0.3004 

143.4292 

0.2524 

25 

7.7124 

0.2522 

0.0002 

0.5044 

0.3084 

179.2865 

0.2524 
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SEC -NO. 

LIQ.T. 
STATIC 
( DEG.R) 

LIQ.P. 

TOTAL 

PSIA 

LIQ. ENTHALPY 
STATIC TOTAL 

8TU/L8 BTU/LB 

LIQ. DENSITY 
STATIC TDTAL 

LBS/Cll. IN 

RE GAS 
FILM 

PR GAS 
FILM 

RE LIQ 
BULK 

RE LIQ 
film 

PR LIQ 
FILM 

1 

56.0 

1152. 

-14.52 

-14.44 

0.00250 

0.00250 

0.983E 

06 

0.662 

0.835E 

06 

0.652E 

05 

0.666 

2 

59.5 

1152. 

-4.39 

-4.26 

0.00243 

0.00243 

0.113E 

07 

0.664 

0.101E 

07 

0. 744E 

05 

0.666 

3 

63.2 

1151. 

6.64 

6.83 

0.00235 

0.00235 

0.130E 

07 

0.666 

0.122E 

07 

0. 822E 

05 

0.665 

4 

67.1 

1151. 

18.73 

19.05 

0.00226 

0.00227 

0.153E 

07 

0.669 

0.151E 

07 

0.932E 

05 

0.664 

5 

71.2 

1149. 

32.08 

32.65 

0.00217 

0.00217 

0.182E 

07 

0.672 

0.190E 

07 

0.1 08E 

06 

0.664 

6 

75.4 

1147. 

46. 83 

47.95 

0.00206 

0.00206 

0.220E 

07 

0.675 

0.245E 

07 

0.126E 

06 

0.663 

7 

77.5 

1145. 

54.61 

56.27 

0.00199 

0.00200 

0.246E 

07 

0.676 

0.284E 

07 

0.140E 

06 

0.663 

8 

79.7 

1141. 

62.69 

65.18 

0.00193 

0.00194 

0.273E 

07 

0.679 

0.330E 

07 

0.149E 

06 

0.662 

9 

81.7 

1135. 

70.76 

74.73 

0.00185 

0.00187 

0.307E 

07 

0.681 

0.390E 

07 

0.161E 

06 

0.661 

10 

82.5 

1131. 

74.56 

79.75 

0.00181 

0.00183 

0.327E 

07 

0.682 

0.430E 

07 

0.1 67E 

06 

0.669 

11 

83.2 

1127. 

77.48 

83.95 

0.00177 

0.00180 

0.342E 

07 

0.682 

0.464E 

07 

0.1 66E 

06 

0.659 

12 

83.6 

1122. 

79.66 

87.45 

0.00174 

0.00178 

0.355E 

07 

0.684 

0.497E 

07 

0.163E 

06 

0.657 

13 

84.1 

1114. 

82.83 

92.69 

0.00169 

0.00174 

0.368E 

07 

0.686 

0.543E 

07 

0.156E 

06 

0.655 

14 

84.5 

1108. 

84.98 

96.07 

0.00166 

0.00171 

0.373E 

07 

0.687 

0.570E 

07 

0.1 56E 

06 

0.654 

15 

85.1 

1101. 

88.05 

99.52 

0.00163 

0.00169 

0.373E 

07 

0.687 

0.583E 

07 

0.153E 

06 

0.653 

16 

86.4 

1090. 

93.93 

105.53 

0.00158 

0.00164 

0.365E 

07 

0.687 

0.595E 

07 

0.162E 

06 

0.655 

17 

88.8 

1077. 

103.56 

113.62 

0.00153 

0.00158 

0.331E 

07 

0.687 

0.580E 

07 

0.149E 

06 

0.653 

18 

92.8 

1063. 

119.21 

125.89 

0.00147 

0.00150 

0.279E 

07 

0.687 

0.520E 

07 

0.140E 

06 

0.652 

19 

100.6 

1052. 

149.37 

151.75 

0.00133 

0.00135 

0.207E 

07 

0.687 

0.395E 

07 

0.122E 

06 

0.654 

20 

109.1 

1048. 

182.06 

182.94 

0.00119 

0.00120 

0.158E 

07 

0.686 

0.299E 

07 

O.UOE 

06 

0.657 

21 

115.0 

1047. 

206.15 

206.64 

0.00110 

0.00110 

0.135E 

07 

0.685 

0.252E 

07 

O.IOIE 

06 

0.658 

22 

120.1 

1046. 

226. 18 

226.63 

0.00103 

0.00103 

0. 131 E 

07 

0.685 

0.241E 

07 

0. 1 05E 

06 

0.659 

23 

125.8 

1045. 

247.87 

248.39 

0.00096 

0.00096 

0.133E 

07 

0.685 

0.242E 

07 

0.U9E 

06 

0.660 

24 

133.2 

1044. 

274. 84 

275.45 

0. 0008 8 

0.00089 

0.135E 

07 

0.684 

0.241E 

07 

0.137E 

06 

0.662 

25 

143.1 

1043. 

309.66 

310.39 

0.00080 

0.00081 

0.138E 

07 

0.683 

0.237E 

07 

0.161E 

06 

0.665 


SEC. NO. R OVER T TUBE DELTA 



TUBE WALL 

STRESSES 


AVG. METAL 

YIELD 

COEF.OF 

YOUNG'S 


LONG. A V 

LONG. OUT 

LONG. IN 

TANG. OUT 

TANG. IN 

TEMP. 

STRENGTH 

L1N.EXP 

MODULUS 

(IN.) 

(PSI) 

(PSI) 

(PSI ) 

( PSI) 

( PSI) 

DEG.R 

(PSI) 

(IN/IN) 

(PSD 


1 

17.321 

0.00732 

-273851. 

-284587. 

-263114. 

59315. 

-21223. 

1436. 

19633. 

0.942E- 05 

0.215E 

08 

2 

15.707 

0.00657 

-273220. 

-286084. 

-260357. 

59890. 

-25903. 

1431 . 

19707. 

0. 942E-05 

0.215E 

08 

3 

14.094 

0. 00600 

-276385. 

-291718. 

-261051. 

61932. 

-32218. 

1464. 

19234. 

0. 945E-05 

0.212E 

08 

4 

12.480 

0.00539 

-278821. 

-297555. 

-260087. 

64453. 

-39261. 

1491. 

18835. 

0. 947E-05 

0.210E 

08 

5 

10. 866 

0.00477 

-281402. 

-304927. 

-257877. 

67844. 

-47563. 

1523. 

18379. 

0.9496-05 

0.207E 

08 

6 

9.256 

0.00417 

-284791. 

-314841. 

-254742. 

73070. 

-58390. 

1567. 

17717. 

0.9526-05 

0.2036 

08 

7 

8.446 

0.00385 

-286217. 

-320967. 

-2 51466 ■ 

76058. 

-64598. 

1588. 

17406. 

0.954E-05 

0.200E 

08 

8 

7.639 

0.00360 

-289577. 

-329011. 

-250144. 

81451. 

-73752. 

1642. 

16587. 

0.957E-05 

0. 1956 

08 

9 

6.833 

0.00334 

-292348. 

-337811. 

-246884. 

87639. 

-84439. 

1698. 

15721. 

0.9616-05 

0.189E 

08 

10 

6.429 

0.00322 

-293877. 

-342824. 

-244930. 

91402. 

-90935. 

1738. 

15089. 

0.964E-05 

0.1856 

08 

11 

6. 107 

0.00317 

-295302. 

-346331. 

-244272. 

95835. 

-97862. 

1796. 

14162. 

0.9686-05 

0. 1796 

08 

12 

5.848 

0.00316 

-295679. 

-347840. 

-243517. 

99968. 

-104288. 

1863. 

13032. 

0.9726- 05 

0.1726 

08 

13 

5.593 

0.00320 

-293013. 

-345649. 

-240377. 

103391. 

-110294. 

1962. 

11305. 

0.9786-05 

0.1606 

08 

14 

5.549 

0.00322 

-291591. 

-344917. 

-238264. 

102445. 

-110326. 

1987. 

10844. 

0.9806-05 

0.1576 

08 

15 

5.593 

0.00329 

-289706. 

-342433. 

-236978. 

100567. 

-108559. 

2013. 

10369. 

0.9816-05 

0.1546 

08 

16 

5.869 

0.00333 

-293152. 

-347717. 

-2385B6. 

93558. 

-100466. 

1955. 

11429. 

0.9786-05 

0.1616 

OB 

17 

6.562 

0.00388 

-288292. 

-338302. 

-2382B2. 

81036. 

-85316. 

2025. 

10141. 

0.9826-05 

0.1526 

08 

18 

8.107 

0.00472 

-289334. 

-333764. 

-244904. 

64166. 

-63515. 

2004. 

10530. 

0.981E- 05 

0. 1556 

08 

19 

11.921 

0.00642 

-292652. 

-3247 89. 

-260515. 

47028. 

-38479. 

1874. 

12842. 

0.9746-05 

0.1716 

08 

20 

16.920 

0.00822 

-287748. 

-310297. 

-265199. 

37446. 

-21965. 

1718. 

15406. 

0. 9646 - 05 

0.1876 

08 

21 

20.441 

0.00964 

-283927. 

-3 01594. 

-266260 . 

35228. 

-15531. 

1674. 

16107. 

0.9626-05 

0.1926 

08 

22 

21.517 

0.00984 

-280547. 

-297727. 

-263368. 

33727. 

-12899. 

1632. 

1674 7. 

0.9596-05 

0.1 966 

08 

23 

21.517 

0.00940 

-275641. 

-294053. 

-257230. 

31788. 

-10982. 

1575. 

17596. 

0.9566-05 

0.2026 

08 

24 

21.517 

0.00890 

-268949. 

-288860. 

-249039. 

29329. 

-8549. 

1511. 

18543. 

0.9526-05 

0.208E 

08 

25 

21.517 

0.00833 

-259808. 

-281557. 

-238060. 

26181. 

-5437. 

1438. 

19602. 

0.9486-05 

0.215E 

08 


TOTAL COOLANT HEAT PICKUP = 5413.36 BTU/SEC 

TOTAL NOZZLE WALL AREA = 1344.24 SO. IN. 

TOTAL NOZZLE WEIGHT « 254.02 LBS 

TOTAL COOLANT TEMP .DELTA = 87.40 DEG. R 

TOTAL COOLANT PRESS. DELTA = 116.72 PSI 

TITLE 

RESULTS 


FUTURE REVISIONS 
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DESIGN PROGRAM 


c ..... 

L DESIGN PROGRAM FOR CDNVbCT 1 VELY LOJLED N3Z/ LtS 

C ..... 

COMMON UL0M,QLIGS,TGE, TGF, TWGG,TWL, TL,PL, VLSF,EML,AL,DF,DM,FF,HGF, 

1 HLF,Lll,C22,C33,QL,ON,RO,RT,HCP, WCp, UHL, WALT, TOTriT , TLS,PT, 

2 HLS,HLT,RS2,KTL,RG»PG,RLl),RL»PF,RUVT»UELT»51GLON»SIGIUO» 

3 SlGLDl»SIGOUT,SlGIN,TMET,SIGALM, ALPHAl.E 

A / INPUT/ AX t 50 > t OHG( 50 ) , AN WE B l 50 ) , PH 1 < 50 I , TW( 5C ) , TWG ( 6U ) , 

5 TGS( 50) , PGS(50) »COEFG( 50) .CUEFLl 50 > , RC( 50) , QRAD* 50 > t 

6 QgEN ( 5C ) , SHEL TH( 5'j) 

CUMMON/ INPUT/ DEN, SEC , G f P I , PGS I » CPH f AMU, ES , WGHT , N , R 1 
CUMMJN/ STATtl/SSI 50 ) / ST ATE2 /Ul 3 ) 

DIMENS1UN LINE* 2) .LIST* 60) .PRINT* 60) 

DIMENSION SEt( 50» 39) • D/150,20), NAME ( 2C 1 . CIKD), OUT<5) 

DIMENSION Cl (01 
EQUIVALENCE (C.QCOM) 

EQUIVALENCE 1RF.RHULF) 

EQUIVALENCE t DPMOM t l)M ) * IDF.DPFR1, { DP » DP TUT ) » * R 3 , RHQLd 1 

EQUIVALENCE (MCAL.Hl 1G> I 

EQUIVALENCE (RL.RCLF) »{RG»REGF) , ( PF.PRLF) , <PG,PRGF) 

EQUIVALENCE (RS2.RHDLS) 

EQUIVALENCE ( PL T , P T ) , I TL , TS 1 
EQUIVALENCE { 1)2 , A X ) 

LUMMDN/ST/H( 10) /MET /TEMPI, TEMP2.X* iC.A) ,XZ(5D> 

EQUI VALL.mCE ( H ( 2 I .CPI . (H( 3) .AM) , ( H I A 1 , AK I , ( il( 5 1 , SVVL I , ( HI 7 1 , R 1 


NAME LliT/CUNSj/CTtX, TEMPI, TEMP2, FI ,AKR, RJUGH, ARLACR, Cl, C2,C3,YM 
NAMEl I sT/iTUFF/ Tw, DUG, AX, TGS.PGS, AN WEB , T WG , COE FI , C JLFG. Rc.PHl , 
lsHtLTri,wGEN,QRAL),CT , WG , WL , PLIN, f L I N , AM J , .JR AZ , OF \ , XT H, YM , I UT , 
2AKR, TEMPI, TFMP2, N TRIAL, ARE ACR, Cl, C2,C3, LIST, L TUP, Fl .ROUGH, 

)X,U IAM.ES, NUMB, COMP 

DATA CT/3. jE- 6, 0.?CP2, I . J, ?. 1 ,0.301 ,3. OO: 1 , D^l ,0.1/ 


REAL KMET 


C 1 = KME T 2 = SIG» 3 = EL, A= APLHA 

G2*2.*i2.17A*12. 

P I “ 3 • 1A1:>927 
G = * 5 

R AU=5 7.29578 
AG«G*G* 

XJ=77b. 161*12. 

LT0P=6« 

DU 1 K = 1 , 60 

1 LISTUI=K 

UO 2 R = 1 » 5 7R 

2 AX(K)=o.o 
KTrt=0 
NtR 1 AL=2u 

A READ { 5 , jT UP F ) 

1 L=TL in 
PL-PL 1 N 
U( 2 ) =CUMP 
NXX=NUMrf*l 

T LhNG= AbS ( AX(NXX)-AX(21 > 

ARLU.\ = U<-IG{ NXX) **2/UlAM*«2 
ARu I V = UHG I 2 ) **2/i) I AM* *2 
WRITE (6,6) 

6 FORMAT Urll I 

CALL CUOLtO.C, 0.9,0) 

WRITt (6, CONST) 

NPL=u 
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WRITE 16,8) 

8 FORMAT (lHl///lX,8M:>tC*NQ. ,3X,12HAxIAL L ENGTH f A X , 9HE NG . Ul AM. , 5X , 

1 6IINU. TUBES, 3X,9HANGLE P H I , AX , 1 1 H TUBE THICK .»6X,7HT.MET.G» 

2 6X, 1JH6TAT. GAS TEMP , AX, 1 AHSTAT . GAS PRfc6S//16X»5H( IN.)«9X, 

35H1 IN. ),20X,6H(l>EG.), 8X,5H1 I N . ) » 9X, BH1 UEG. Kl , ilX,8H(DEG. R) ,9X, 
46H1PSIA)///) 

DO 1C J=2 » NXX 
JJ=J-1 

ID WRITE (6,12) J J , ( DZ( J , I i , 1=1,0) 

12 FORMAT (2X, 13 , FI 6. A , FI 5. A, F12. 1 , F13. 3, F14. A, F15. 3, F19. 1 , F16. 2) 

WRiTt 1 6 , 1 A ) 

IA FORMAT (1H1///1X»9H SEC. NO. ,2X,9HGAS C3EF . t 2X , 1 DHCOGL . CUtF . , 1 X, 

11 IHKAO .OF CUR V, 3X , 6HU RAO. , 6X , 8HNUL .HEAT, 5X, 12HSHELL THICK.// 
237X,t»H{ 1 4 • ),6X, 7HQTU PER,5X,7HOTU PER , 9X, 5H( I N . ) / 

348X,9H>U. IN-$EC,3X,9HCJ.1N-SEC///) 

00 16 J=2,NXX 
JJ=J-1 

16 WRITl (6,18) JJ,(OZ(J,I), 1=9, IA) 

18 FORMAT (3X, l 3 , 2X , 2F1 1 . A , 61 A. A , F 11 . 3 , 2F13.3) 

WRITE (6,20) OEN, 6S, AM J, BRAZ 

20 FORMAT 1///1X, 18HAVG.ME1 AC OENS l T Y= , F6. J , 2X , 9HLBS/CU. I N, 3X, 

1 22HSHELL YOUNG'S MO0‘JL JS = , E 1 0 . 3 , 2 X, 3HPS I , 3X , 16HP0 I SSU i* S RAT 

210=,i-6.3,3X,12Hi3RA2E THIGK=»F7.4»2X»5H( IN.) ) 

WRilt 16,22) WG,WL,TGMNXX1 , TL,PGS(NXX) , PL , 0 1 AM, TLENG, ARCON, 

LARUlV 

22 FORM A ! I1H1//9X,21HPR0PlLLANT FLOW RATE ,F9. 4, 2X, 7HL0S/SEL, 23X, 

I 17HCU0L ANT FLOW RAT t , 4 X , F 9. 4 , 2X , 7HL 8 S/ SEC//9X , 1 9HCHAMBE R TEM 

2PEKA1 URt, 3X, F8. 1,2X, 6HJEG.R , 2AX, 22HEN TRANCE COOLANT TEMP., 
3F8.2t2A,bhUfcG.R // )X, l 6HCHAMBER PRESSU T E , 6X, F8 . 1 , 2X, 6HPS I A , 
424X,22HtNTRANCE COOLANT PRE SS , F 3 . 1 , 2X , 6HP S I A ///A3X, 15HTHKUAT DIA 
6METER, iiX»F9.4,2X»5H( IN. )/A3X, 1 9HT U TAL NOZZLE LENGTH, 8X, F 10. A,2X, 
65H1 IN. )/A3X,26HCDNVERGtNT AREA RATIO A/ A* , 3X, F8. A , /A3X , 25H0I VERGEN 
71 AREA RATIO A/ A * , AX , F8 . A/ / ) 

M T = NT RIAL 
DM=C .0 
UF=C . 0 
UP=0.0 
0TUT=0.0 
A TU T = 0 .0 
WGrtT = 0 •(. 

C PROGRAM MAIN LOUP 

00 128 N= 1 , NUMB 
NSEX=NUMB 
00 2 A L=l, 14 
J=N+KTH+i 

2A UZ(1,L)=0Z(J,L) 

CPH=LU5(PHI/RA0) 

Rt=G*U»1G 
R Z = 0 . 0 
NPL = 3 
KA=G 
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c 


CALC. OF GAS SIDE HEAT FLUX 


26 DU 2 8 M= l * MT 

T Gt=TGS+AKR» ( TGT-TGSI 
T GF=G* rWG«-.28»TGS+.22*TGt 
F=2.*TW+t*RAZ 
OG1- JG2 

WG2=wGtN»SrttLTH 
UN- (Jo 2 

QZ=G*(wGl+jG2» 

CALL LOUL (TGS,PGS,4» 

RHUGi=R 

CALL COUL ( TGF, 203,4 ) 

RHOgF=R 

RG=WG/.7B5/OrtG/AM*RHOGF/RHJGS 

AMu=AM 

NCT = 1 

R=RG* » . 8 

AKG=AK 

CPG=lP 

PG=CF*AM/AK 

AKR X = AKR 

AKR=PG**. 33333 

IF { Aui( AKRX-AKR) .LE.CT (bll GO TO 30 
28 COlMl 1 NUfc 
30 p =PG» * • 3 

rtGF=CUbFG«AK»R*P/OHG 
UC=HGF»{ TGE-TwGI 
UCUM= UC+U RAO 

C CALC. OF COOLANT SIDE WALL TEMPbRATURt 

I MET = T wG 
P 2 = P I / AiSIWEti 

RT=(P2*(RQ+TW*CPHl-F/2. )/< l.-P2*CPH) 

RT0 = R1 + T W 

RAjlrt = KlU*ALOG< RTO/RT ) 

A=OCuM*R AUT w 
JO 3g K= 1 , MT 
lK=Khtl(l, IKfcT) 

TWL=T WG-A/T* 

IF (TWL-TL) 32, 32,34 
32 TWG=TwG*100. 

GO 10 26 
34 T =T ML r 

TMLT = G*( TWG«-TWL ) 

IF(AUj(T-TM£ 1 ).LT.LT<3) ) Gj 10 33 
36 CONTINUE 
38 C2^=i. 

WL I - WL/ AN WE 6 
RZ=(RT/RC) **2 
R l=RL + CPHM RT + TW) 

WCP = 2 . *R T 

K 0= 1 

H4=PI*KI*RT 
1 F ( N. G i • 1 ) GO TO 43 
HLP=wGP 
GO TJ 42 
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CALC. OF COOLANT ENTHALPY INCREASE 


L 

AD 00= ( Ul-DHG 1/2. 

0X=AX1-AX 

SfcC=SUKT( 00*00+ OX *DX> 

WALL = UU»ANwFBl/A>jWI:e + Rl+TWl + Trf)*PI*SEC*G*AREACR 
walt=wall*anwfb 

ASd=G»(ANWfc31/ANWE0*UHG/01+l.)*SfcC*WCP 

UG=02*A.>rf/WLT 

U A=G* I uCOM+Ol )*WALL/WLT+QG 

OL=OA*WLT*ANWErt 
HLT = HL Tl + OA 


C COOLANT SIOE CALCULATIONS 

A2 NCT=1 

00 112 NH= i » MT 
K91 = j 

IFfN.til.il GO TO AA 
CALL CUOL I TL»PL» A) 

R6=R 
AB=AM/R 
GO TO A6 

AA CALL CUUL (HLl*PLT»-3) 

TL=SS( 7) 

IFINH.EO.IJ HLS=HLil+QA 
CALL CUUL ( HL S » PL » - 3 ) 

TLS=SS ( 71 

CALL CUOLI TLS,PL»A) 

RB = R 
A ri = AM/ K 

RLB=HLT/AL*OUL/AM 
A6 CALL COOL! TrfL.PL, A) 

AW=AM/K 

HLF = wCOM*RTO/RT/( TrfL-TL 1 
JL 1 (J.* = uCUM*R 1 0/ RT 
TF=TLF 

TLF=b*(TwL+TL) 

CALL CUUL ( TLF, PL » A ) 
kF=r 

AML =ttM- 

CPL=LP 

AKL =AK 

PF=CP*AM/AK 

S = *lLT/AM*R/RB 

Cll=l.+.01A57*Aw/AU*Cl 

C22=l. 

IF(N*NM.cQ.l) GtJ TU A8 
RL=wlL 1 / AL*OHL/AM*RriOLF/ftHULB 

1 FI RC*C2 • bU. f- • D) GO TU AB 
IF(RL*K2.LT.6.^) GO TO AB 
1FIRC.LT. D.: ) C22=l . 

IFIRC.Ol.O.C) C22=(RL*R21*»0.05 

A8 CC=C 1 1 *C22 

CSC=tUtFL*AK/HLF*S**J.6*CC*PF**D.A*1.5U 
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C GEOMETRY CALCULATIONS 

DO 54 J=1,NTRIAL 
HI = 2 • *R T FWCP 
H2=CSu/Hl 
H3-P 1 *R T + HCP 
H5=H4/Hl 

UU t>0 K= 1 1 NTR IAL 
HX=HCP 

WP=H3+i.*HLP 
HCP=H2«rfP*»3.2-H5 
IF(HCp.Lfc.O.O) GO TO 9 A 
I F ( Aus(HX-HCP) ,LT.CT(2M GO TO 52 

53 CONTINUE 

52 RU=R1+CPH*HCP 
rfX=*CP 

NCP=2.*P2*RD-F 

lF(AbS(rfX-WCPl.Ll.CT(2) ) GO TO 56 

54 CONTINUE 

56 AL=G» ( H4*HCP*H1 ) 
uHL=4.*AL/WP 
lFlN.tu.il GU TO 114 


C CALC. OF STATIC PRESSURE DROP 

IF (KA.EO.O) GO TO 06 
58 A V=G* ( AL 1 + AL ) 
irfPV=G* ( rfPl+WP I 
0V=4.*AV/rfPV 
GV= rfL T / AV 
TV=G*t TLS + TLSl) 

P V=G» { PL+PL 1 ) 

CALL CJUL( TV, PV, 4) 

R=GV*UV/AM 
C33 = l . 

IF(KL*C3.EU.0.3) GU TO 63 
I F { R*R 2 . L T .6.0) GU TO 60 
C33=(K*R2) **0.^5 

63 F2 = F 1 

1F( RuuoH. to. 0.0 ) GU TU 62 

FF = ( l.G/I-4.0*AL3G10lRUUGH/l3V/3.7+i.253/R/FL**C.5 1 >1**2 
GU TO 64 

62 FF= Il.u/ (4. :*ALJCI3(Fl»*j.5*R»-3. 40 »)**2 

64 F l = Fh 

1F ( Aus( F2-FFI .GT.Crm ) GO TO 63 
FF=FF*L33 

r =2 .*ov*gv*sec*ff/ag /ov 

A1=I./ALI/RS1 
A2= 1 . / AL 
A=WLT*hLT/AC/AV 

00 74 J= I , MT 
UU 70 K= i , MT 
CALL lOUL (HLS.PL, -3) 

RS2 = R 

VLS=wL1/KS2/AL 
R V = G* ( KS 1 *RS 2 ) 

OM=A»( a2/RS2-AI ) 

UF=T/KV 
UP = DF + DM 
PLX=PL 
PL = PL 1-L)P 

IF(PL.GT..3*SFEl 1,8) ) i»U TO 68 
PL=PL 1 

T WG= T Wo 1 -2 33 . 

NPL=NPL+1 

IF (NPL.bT.il GO TU 66 
GO TU 26 
66 TWG=T*>lt>l 

IF (NPL.GT.2) GO TU 13, 

GO TU 26 

60 IF (AiibtPLX-PLl .LE.CTI4M GO TU 172,1181, NOT 
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C EXIT POINT FDR FINAL CONVERGENCE 

iNCT = 1 
PLX=PL 
70 CLUTlNUE 
72 HLX=HL b 

HLS=HLT-VLi»*2/G2/XJ 
HTUL = Ail:>(HLX*CT( 6) ) 

IFI AbiIHLX-HLS) .LE.HTUL > GO TO 76 
74 CONTINUE 

C CALC. OF TOTAL PRESSURE ANJ TEMPERATURE 

76 CALL COOL (HLS,PL»-3> 

TLS=SS(7) 

SVL=SVVL 

EML-VLS/SVL 

TV=G*(TLS+TLS1) 

PV=G*l PL+PL1 ) 

CALL CuOL(TV,PV,4) 

RSA=R 

UO 78 L=l,MT 
TTV=G*( 1L+TL1) 

PLTV=G*(PLT*'PLTI) 

CALL COOL ( TTV,PLTV,4) 

RTA=R 

OPT =XJ *KTA *^A*{ L.-T TV/TV l -0F»RTA/R6A*( T TV/TV) 

PLTX=PLT 
PLT=PLT1+0PT 
T LX= TL 

CALL COOL ( HL T » P L T » - 3 ) 

RTL=R 
TL=SS ( 7 I 

IF( A6SIPLTX-PLT J .LI.CTt 4>.AN0.ARS( TLX-TL) .Lt.CTIS)) 30 TJ 83 
78 CONTINUE 

U MACH NO. LIMITATION PROCEOURC 

8C IF ItML.LT.YM) GO TO ( 92 . 82 , 96 ) , KD 
82 IF (XU.LT.3I KO=KD+l 
VLS a YM*3VL 

HLi a tlLT”VLS**2/G2/XJ 
CALL COUL(HLStPL,-J» 

R S 2 = R 

AL = WL I/Ri2/VLS 
Z=WLr/AL*OHL/RHOLa 
00 8*» K-l » 6 
rtCP a 2./Hl*(AL-H4*G) 

RO=R I +CPH*HLP 
dX=dCP 

WCP-2.*P2*R0”F 

WP=PI*KT+2.*HCP+hCP 

JHL*4.*AL/rfP 

lF(ABi(HX-NCPI.LT.tT<2) ) GJ TO 96 
84 H1=2.*KT+WCP 
GO TU 9& 

C AREA MODIFICATION FOR TUBE SPLICE 

86 IF ( ANrfEd-ANWCdl J 89,58,88 
88 K A- l 

NR=ANrttBl/ANWEd 
ALi*wR*ALl 
RT1 = WR*KT 1 
«CPl=WK*rtCPl 
rIll*^.*RTI+vJCPL 
H4I«P I »RT 1**2 
00 9- K = l,5 

HCPl=2./HlL«( AL1-H41*G) 

RU1 -R I I+uPHl*HCP I 
WXL=WtPl 

rtCPL-2.*PI/ANvlFb*R01-( 2 .*T4U-BRAZ) 

NPI-Pi*RTl+2. *HCP I+WCPl 
IF (rtdj(WXl-wCPI) .LT.CI (21 ) GO TO 68 
90 HI W.»RT1 + WCPL 
GO TO SB 

92 IFUu.uT.RlJ GO TU 112 
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c 


AREA MODIFICATION FOR MINIMUM AREA 


9A AL = H<t/2. 

RU = RI 

DHL = 1 . 222*RT 

HCP=C.O 

UCP=2.*RT 

C EVALUATION PROCEDURE 

96 Z*WLT/ AL * DHL /R MO LB 
DO IlU K91 = 1 » MT 
UAI =UA 
TF=TLF 

CALL COUL (TLF.PLiA) 

RHULF=R 
RL=Z/AM*R 
PF=CP *AM /AK 
CPL-CP 
AKL=AK 

CALL CUOLl TLS.PL.Al 
RB=R 

RLtt=rfLT/AL»OHL/AM 

ab=am/k 

CALL COOL! THL.PL, A) 

AW=AM/R 

CU = 1. + .01A5 7*AW/AB*C1 
CC=Cli 

IF(RC*C2.EU.O.OI GO TO 96 
IF(RL*RZ.LT.6.')> GO TO 98 
C22=(RL*RZ) »*0.05 
IF(RL.LT.O.O) C 22 = 1 • 

CC=CC*C2 2 

98 HLF=CUtFL»AKL*RL»» j. 8*PF**3. 4*CC/DHL 
TWL=UCOM»RTO/RT/HLF+TL 
TLF=G» ( TWL+TL ) 

DO 106 K= l ? MT 

A=UCOM*RADTrf 

T=TMEI 

TMET=G*I TWL+TWG) 

THG = A/KMtTU,TMET)+TWL 
TGF=G*TWG+ .28*TGS+.22*TGE 
CALL CUOL ITGF,PGS,4) 

RHUGF = R 

RG=HG/ . 766/ DHG/ AM*RHUGF/RHJGS 

AMG= AM 

AKG=AK 

PG-CP *AM /AK 

HGF = CUfcFG*AK*RG*O.8*PC*»0. 3/DHG 
UC= HGF* l TGE - TWG ) 

07=QL0M 

dCOM=uC+URAD 

UT0L=0C0M*CT(7> 

IF IAtt^UCDM-07) .LE.QTDL) GO TO 103 
GO TO 106 

103 IF { Ad*l TF-TLF) .LE.CTI3M GD TO 102 
GO Tj 110 

102 IF ( Abil I-TMETI .Lt.CTI 3) ) GO TO 10A 
GO TO i06 
iG* NCT =2 

JL lQb=wCUM«RTO/ RT 
GO TU 56 
106 CONTINUE 
GO TU 58 

108 UC0M=G*(UC0M+Q7) 

UA=G*(gLJM+Ul > *WALL/*LT*UG 

UL = vjA*RLT*ANWE8 

IFIN.hU. 1) GO TU 113 

HLT = ML T 1 +QA 

CALL COOL (HLT,PLT,-3I 

TL=S5( 7) 

110 CONTINUE 
112 NCT=2 

GO TO 116 
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C CALC. DF TOTAL PRESSURE AT INLET 

114 CALL COOL l TL , PL , 3 I 
VL=(rfLT/R/AL)**2/G2 
PT=PL+VL*R 
CALL CUOL { TL , P T , 3 ) 

HL T-H ( b ) 

HLS=HLT _ VL/XJ 
DO 116 K= 1 » MT 
CALL COOL (HLS, PL, -3) 

VLS=WLT/R/AL 

RS2=R 

PT=PL+VLS**2*R/G2 
CALL COOL <TL.PT, 3) 

HLT=h( o) 

HLS=HL T-VL S* *2/ G2/ XJ 
CALL COOL < HLS , PL , -3 ) 

TLS=Sbl 7) 

SVL=SVVL 

EML=VLS/SVL 

I F( ABS l HX-HLS I •LT.CTC6)) GO TO 118 
HX = HL S 

116 CONTI NUfc 
GO TU 118 

118 CALL STRESS 

C END OF PROGRAM CALCULATIONS 

C RESULTS REUUIRED FOR NEXT STATION CALCS. 

PGS i=PGS 
A L 1 =AL 
01 = UHG 

AX 1 = AX 
UHL 1 = JhL 
ANrtfcb 1 =AN WEb 
RT 1 =RT 
TR1=TR 
HLT 1=HLT 
T LS 1 = TLS 
PLT 1=PLT 
WCPI=WCP 
R I l=R I 
CPH1=CPH 
PH 1 1 = PH 1 
RP1=UP 
TL1=TL 
T RG1 = T wG 
TrfGG=TWG 
TWGl N*1 ) =TWG 
PL l -PL 
AKT = 1 MET 
U1=QCUM 
RS=RS2 
HLS 1 = HL S 
RSI =R S2 
T UT HI =RT +HCP 
jTUT=gTOT+UL 
ATOT=ATUT+WALT 
VLsF=VLS/12.9 
00 12u K= 12, 5P 
JJ=K-11 

1 2 J StE< N» J J I =Ct K } 

AKL=G.C 
PF=0 .0 
RL= J.O 
AMG=J.O 
AKG= J. 0 
PG=0. 

RG=0 . 0 
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OUTPUT 

IFIN.NE.l) GO TO 124 
WRITE (6,122) 

122 FORMAT ( 2X.7HSEC. NO. ,2X,7H0/A 0uT,4X,6HU/A IN, 3X, 9HEFF.GAS T,2X, 
16HT.REF.,4X,7HT.MET.G,3X,7HT.MET.L, 5X,6HL1U.T., 

24X,6HL 10. P. ,4X, 7HL I Q. VEL , 4X, 7HM AC H N0,3X,9HTUBE AREA/ 
374X,5HT0TAL,5X, 6HSTATIC/ 

410X,18HltJTU./SQ. IN. -SEC. ) ,4X,7H(DEG.R) , 3X, 7H( OEG. R ) , 3X, 7H( OEG.R ) , 
53X, 7H( UEG.R) ,5X, 7HI OEG.R) ,3X, 6H(PSIA) ,4X, 6MFT/SEC, 15X, 

68 H( SU • IN*)///) 

124 WRITE (6,126) N , ( C ( K ) , K* 1 , 1 1 ) 

126 FORMAT (3X,I3,2F11.2,4F10.0,F11.1,2F10.3,F12.3,F13.4) 

128 CONTINUE 
N=NUMB 
GO TU 132 

133 N=N-1 

WRITE (6,131) 

131 FURMAT (///25X,31H*» PRESSURE BELOW MINIMUM VALUE) 

132 WRITE (6,134) 

134 FORMAT ( IH1///2X, 7HSEC • MO . , 2 X, 10HFR I .P .OR OP, 2X, 1 OHMOM . P . DROP , 2 X, 

1 6HFRI.FACT, 7X.3HHGF, 1 IX, 3HHLF, lOX, 2HC1, 9X, 2HC2 , 8 X , 2HC3, 7X , 

2 9H(J/SECTI0N,6X,7HQ/A NUC// 1 3X, 5H( PS I ) , 8X, 5H ( PS I ) , 1 7X, 2 IH( 6TU/SU.IN 
3.SEC.0EG.R), 39 X, 7HBTU/SEC, 4X , 13HBTU/SQ. IN. SEC///) 

00 136 J = 1 , N 

' 136 WRITE (6,138) J,(SEE(J,I), 1=1,13) 

138 FORMAT ( 3X, 1 3 , 2F 1 2 . 2 , 3F 1 3 . 5 , 3F 1 1 . 3, F 12. 1 , F 1 5 . 3) 

WRITE (6,140) 

140 FURMAT( 1H1///2X, 7HSEC.M0. ,6X, 9HSHELL RAD,8X, 

1 llHTUaE RADIUS, 6X, 11HTU8F HE I GHT , 6X , 10HTUBE WIDTH, 9X, 7HHYD.D 

2 I A , 7X , 9HWALL AR E A , 6X , 1 7HT0 T A L TUBE HE I GHT // 1 8 X , 5H ( I N . ) , 1 1 X, 6H( I N. ) 
3, 12X,5H( IN. ) , 13X,5H( IN. ) ,12X, 5H( IN. ),1DX, 7H1S0. IM) , 12X,5M( IN. )///) 

00 142 J = 1 » M 

142 WRITE (6,144) J,(SEE(J,I), 1=11,17) 

144 FORMAT (3X, I3,7F17.4) 

WRITE (6,146) 

146 FORMAT (lHl///2)^» 7HSEC . MO. , 2X, 6HL IQ.T. , 3X, 6HL I U . P. , 6X, 

1 I2HLIQ. ENTHALPY, 11X, IIHLIO. DENSITY, 9X,6HRE GAS,5X,6HPR GAS, 

24X , 6HRE L I U , 7 X , 6HRE L I (J , t»X , 6HPR L I U/ 1 1 X , 6HSTA T I C , 3X , 5HTU I AL , 4X, 6HS 
3TATIC,6X,5HTOTAL,6X,6HSTATIC,5X,5HTOTAL,7X,4HFILM,7X,4HFILM,6X, 
44HtJULK,9X» 4HF ILM, 8X, 4HF 1LM/ 11X, 7H( 0EG.R),2X,4HPSIA,5X,6H3TU/LB,6X, 
56HBTU/LB, 9X, 9HLBS/CU. IM/// ) 

00 148 J = 1 , M 

148 WRITE (6,130) J,(SEE(J,I), 1=18,28) 

153 FORMA! (3X,13,F13.1,F9.J,2F11.2,2F11.5,bl3.3,F9.3,2tl3.3,F9.3) 

WRITE (6,152.) 

152 FORMAT! 1 H 1 // / l X , 7HSEC . MO., IX, 8HR OVER T,2X,10HTUBF DELTA, 22X, 

1 iBHTUBE WALL STRESSES, 15X.9HAVG. METAL, 3X,5HYIELD,5X, 

2 7 HCOEF . Uh , oX , 7H YOUNG * S/ 3 1 X , 7HLUNG . A V, 3 X , 8HL0NG . OUT , 4 X , 7HL0NG. IN, 
34X, 8HT ANG. 0 JT» 3X, 7HT ANG . I N , 4X , 5HTEMP . , 3X, BriSTRE NOTH , 4X , 7HL I N. EXP , 
46 X , 7HM0DULUS/ 21X,5H( IN. ),6X,5H(PSI) ,6X,5H(PSl) ,6X,5H(PSI) ,6X,5H(PS 
51) ,6X, 5H( PSI ) ,5X, 5HUEG.R,4X,5HIPSI ) ,6X,7 H( In/ I N ) , 7X , 5HI PS I ) /// > 

DO 154 J= 1 , M 

154 WRITE (6,156) J,(SEE(J,1), 1=29,39) 

156 FORMAT (2X,13,FID.3,F11.5,5F11.0,F9.0,F11.0,2E13.3) 

OELTT= TL-TL I M 
UELTP=PL1N-PL 

WR1TL (6,158) QTOT, ATOI , WGHT, DELTT,DELTP 
158 FORMA T( 1HI////9X,27HT01AL COOLANT HEAT PICKUP = , F 1 3 . 2 , 2X , 7HBT U/ SEC 
1 / /9X , 2 7HT0T AL NOZZLF WALL AREA = , F 10. 2 , 2X , 6 (SQ. I N. //9X, 

227HTUTAL NUZZLE WEIGHT = , F 1 3 . 2 , 2 X , 3 HL BS// 9X , 2 7HT0T AL I.OOLANT 

3 TEMP. DELTA =,F10.2,2X,5HOEG.R//9X,27HTOTAL COOLANT PRESS. DELTA = 
4 t F10.2,2X,3HPSI ) 

WRITE (6,160) 

160 FORMA T I//50X,5HT1TLE/ / 5 OX , 7HRESULT s//// / 50X , 1 6HF JT URE REVISIONS//) 
GO T 4 
END 



SUtiRuUTlNE COOL CX,Y,J) 

CUMMj.M/ilAltl/3S(50)/STATE2/U( 3 > / STATE3/C ( 2 i 5 ) / ST ATE JK ( 5 'I 
COMMUN/SI/HC 10 ) 

COMMON XZ ( 6) » TL 
tUU i VALE xCE <TL,TS) 

tUU I VALENCE CP , S S C 6 ) ) , C T ,SSC26M,*NF , Sit 5) ) , C CP, SSC 31 ) ) 
EUU I VALENCE ( AK » SS ( 30 ) ) , C VL , SS ( 1 5 ) I , C VG, SSC 16) ) , ( VF , Ss ( 19 ) ) 
EQUIVALENCE (AM f SSC27)), CHF,S$C17)> , C NCAL , H( 13 ) ) 

P = Y 

IFCJ.Nb.CI 60 TO 10 

NCAL= I 

MF=1 

U=-2. 

UC 3) = 1. t-36 

VL=.25 

VG=4. 

SS ( b) = • 5 
VF = . 5 
SSC 8 ) = VF 
CALL STATE S 
RE FURiM 

1C NCAL=NCAL+1 

IFCJ.bl. 0) GO TJ ( 11 i 11, 11, 12) ,J 

SSC 7)=TS 

Sit 12 ) = X 

CALL STA1E (J) 

Hlb)=Ss< 12) 

H(7)=1./SSC8) 

H C 2 ) = s s ( 1 9 ) 
rtC5)=ssC9) 

RETURN 
il SS C 7 ) =X 

CALL STATE U) 

HC o ) =SS( 12) 

Ht7)=i./sS<&) 
rt(3)=Sit n 
RETURN 

12 SSC 26 ) =X 

LALL cTATE CJ) 

Hl2)=Ss< 31) 

H ( 3 ) = S^C 27 ) 

H(^)=SS(30) 

MC S ) =SS ( 9 ) 

H( 6 ) =S->C 1 2 ) 

H(7)-1./SSC19) 
rlt 8 ) =sSC 28 ) 

RETURN 

ENO 
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SUBROUTINE STRESS 

CUMMDN QCOM.OLIQS.TGE.TGF, TWGG.TWL, TL.PL, VLiF, EML, AL,DF,DM,FF,HGF, 

1 HLF,Cll,C?2,C33,QL,0N,R0,RT,HCP,WCP,DHL,HALT f TUTHT, 1LS.PT, 

2 HLS,HLT,RS2,RTL,RG,PG,RLB,RL»PF,ROVT,UfcLT, SIGLON, SIuLUO, 

3 SIGL3I,SIGUUT,SIG1N, TMET.S1GALM, ALPHAl.b 

4 / INPUT/4 X < 53 ),DHG< 53) , A, WEB I 53), PH II 5 3) , T.V< 5„ ) , TWGI 30) , 

5 TGS153) »PGS(53I ,CJEFG{ 50) ,CJEFL( 53) ,RC(53) ,yRAlM5j> , 

6 GGENI 50) , SHLLTHI 50) 

COMMUN/ INPUT/ DEN, SEC.G, PI , PGS1 , CPH , AMU, ES, rfGHT, N, R l 
REAL KMfcT 

IF(N. GT. 1) GO TO 9 

RAU=57. 29578 

\/T = 0,o 

VS=0.0 

AT*0 ■ 0 

AS=0.0 

ENS=m .0 

£NS=G .0 

9 E-K.MLT I 3, TMFT ) 

A = XMC.T ( 4 , TMbT ) 

OS 1 -US 
bNS l =bNS 

US=2.*IRI+HLP/CPH) 

bNS=-G*(PGb*PGSl) *t DS1»0S1— DS*DS> 

ENS = tiW( 4.»DS*CPH) +EN:>l*DSl/OS*COS( ( PH I -PH 1 1 ) /RAD ) 
tNl =PGS*US*G/CPH 

U RX=RT+G* T R 
£6=0.0 
t7 = bo 

IF( SrtELTrt.bO.0.0) GO TO 12 
b6=( fcNS-AMu*£NT )/Es/$MELTH 
E7={t.NT-AMU*£NS)/LS/SHbLTH 
12 Z=RX*(PL-Pob)/E/T W 

ALPHA l =KME T ( 4 , G * ( TMET+TL) > 

R=ALPMAl*( TMET-TL ) 
b 8 = W-AMU*Z 
fc 9 = Z 

l)ELT = 2 . * t t E9-AMU* (b6-E8n*RX-Pl/ANwE8*E7*RI> 

T 1 1 = L * ( 1.06B*TR*UbLT/RX**2/( 1 . -AMU* *2 ) - A* ( T RG-T RL > / z . / ( 1 . - AMU M 
I l2 = b * ( Z-. 2796*U£L T *T W**2/ RX**3/ I i.-AMJ**2) ) 

S I GOU 1=T12+TI1 
SlGIN=T12-Tll 
bl5=E6-ba 
S I GL UN = t *E 1 5 

S 1 GMB T = b * A * ( TWG-TrtL)/(2.*( l.-AMU) ) 

SlGLUU=blGL3N-S IGMuT 
SIGLuI=SlGL3N+SIGM6T 
SlGALM=RMET ( 2. TMbT) 

RUVT=KA/TR 

ATu=aT 

AT = T R* ( 2.*HLP+P I*RT ) * ANRbB 
VT=.5*hbL*{ ATU+AT ) +VT 
iU = 2. *RU+SHbLTH*CPH 
A SU = Ai 

AS=P1*S0*SHELTH 
VS=0.3*SLC*{ASU+AS)+VS 
lFlfM.bv.L) RETURN 
RGHT = uE-N* { VI +VS ) 

RETURN 

LNU 

REAL FjNLTUM KMbT(N.TL) 

L OMMU.N / Mb T/IEMP1, 1 1 MP2 , X I 1 3 , 4 I , XZ ( 30 ) 

JlME.MblU>J 1 bN ( 4 ) 

DATA ( 1 EN ( I) , 1=1, 4}/. 0331, 10C0 • , 1 0 j 00 30 • » • DOC 1 / 

T = I E 

I F ( TEMPI. GT. I ) I = TEMP I 
I F ( TLMP2.LT.T ) T=TLMP2 
T = T/T£M2I 
KMET = X ( 8 » N ) 

K = 8 

DO 1J J= i » 7 
X =K- L 

10 KMET=X(K»N)+T*KMET 
KMtT=KMtT*TEN(N) 

15 RETURN 
£NO 


EVALUATION PROGRAM 


C EVALUATION PROGRAM FOR CONVECTIVELV COOLED NOZZLES 

COMMON X(IQ»5)»T£MP1»TEMP2, TEHP3 

1 / INPUT / A X ( 50) ,Dr!G( 50) ,ANWEB l 50) , PH 1 150) ,TWl 50) , AREAL! 50) , 

2 TGS< 50),PGS(50),COEFG(50),COEFL(50J ,AREACR{50) , ROUGH! 50) , 

3 RC( 50 ) » QRAD! 50 ) , UGEN! 50 ) , SHELTH! 50 ) , TKCER !50),TES1C(50) 

A /0UTPUT/DELPF<50) , DELPHI 50 ) , CFLIQ! 50 > , HGF ( 50 ) , HLF i 50 ) , 

5 C1(50)»C2(50)|C3( 50 ) , QSEC! 50 ) , QNUCPT t 50) , RADSH ( 50) , 

6 RADIUS! 50 J , HEIGHT! 50), WIDTH! 50) .DHL! 50) ,AHALL 150), HIGH! 50), 

7 TLSIAT(50) ,PLTOT 150) » HS( 50) ,HT( 50) , ROS ! 50 ) ,RUT!50) , 

8 REGFI50) .PRGF150) ,RE3ULK{50) ,RELF!5C) ,PRLF( 50) ,ROVT!50) , 

9 DELT!50) ,SIGLON!50) ,SIGLOO! 50 ) , SIGLOI1 50 ) , S I GOUT 1 50 ) , 

A SIGIN! 50 ) ,TMET!50) , S I GAL Ml 50) ,ALPHAl!50) ,E(5C) 

B /STATE1/ STORE! 50 ) / STATE2/UN I TS , COMP , CONV/STATE 3/CS! 215 ) / 

C STATEA/JUNK! 50) 

DIMENSION DUMMY 1(50,18), DUMMYO! 50,39) ,CT(9) 

EQUIVALENCE! DUMMY I, AX) , ! DUMMYO, DELPF ) , 

1 1PST3A, STORE! 6) ) ,( CPSTAA , STORE < 31 )),< TSTATA, STORE! 26 )) , 

2 { RST A3,STORt!6) ),( AMUSTA , STORE < 27 } ),! AKSTAA, STORE! 30 )) , 

3 (TSTA3, STORE! 7) > ,( ENTH33 , STORE (12 ) > , ! RSTFA, STURE! 19) ) 
REAL KMET 

LOGICAL LESTCl, LESTC2, LESTC3 , LESTCA, LESTC5 

NAMELIST /STUFF/TW,DHG,AX,TGS,PGS,ANWEB,AREAL,COEFL,COEFG,PHI,QRAO 

1 , TKCER, T tSTC,SHELTH,RC,ARbACR, ROUGH, QGEN, NUMB, COMP, AMU, 

2 TLIN,PL1N,BRAZ, WL , WG, D I AM , DEN , ES , TEMP 1 , TE MP 2 , TEMP 3 , X , 

3 QEST, TCElt , 1 WGE , TMET1 ,F1 , AKR, CT/CONST/CT, X, TESTC, TEMPI , 

A TEMP2,TEMP3,QfcST,TCER, TWGE »TMETi,Fl,AKR 

DATA P I, P 12, A, A 1,A2, A3, AA.CT/3.1AI 5927, 1.5707963, 0.003333333, 

1 9 337.93 2,0.017A5 3292, 1 . 3868612E-7, 5 . 7870370E-A , 3. OE-6 , C. 0002 , 

2 1.0,0.1,0.001,0.0001,0.001,0.1,0.3001/ 

DO 2 1=1,900 

2 AX<n-0.0 

TEMP3=0* Q 
DO A 1=1,10 
00 A J=1 , 5 
A X! I,J)=0.0 
6 WRITE! 6,8) 

8 FORMAT (1H1) 

READ! 5, 10) LESTCl, LESTC2,LESTC3,LESTCA,LESTC5 
10 FORMAT -( 5L 5 ) 

READ<5, STUFF) 

CQNV= 1 . OE-6 
UNITS=-1 .0 
CALL STATES 
QNMi=Q£ST 
QI=QEST 
P=PL IN 
TL=TL IN 
T I = TL I N 
NXX=NUMG+1 
T GT = TbS ( NXX ) 

TLENG=ABS(AX(NXX)-AX(2) ) 

ARCON=DHGI NXX)**2/UIAM**2 
ARDIV = UHG(2) * *2/0 I AM* *2 
MR I Tfc!6, CONST) 

WRITE ! 6, 12) 

12 FORMAT! 1H1///1X, 8MSfcC.NU. »3X,12HAXIAL LENGTH, AX, 9HENG.UI AM. , 5X, 

1 8HN0. TUBES, 3X,9HANGLE PHI , AX , l IHTU bE TH I CK . , AX , 1 3HAK E A PER. T 

2UB£,5X,13HSTAT.GAS TEMP, AX, 1AHSTAT .GAS PRESS//16X, 5H{ IN. ) , 9X, 

35H! IN. )»20X,6H(D£G.),8X,5H( IN. ) , 9X, 8H( SQ. IN. ) , 11X,8H(0EG. R) ,9X, 
A6H! PS l A)/// ) 

DO 1A J = 2 » NXX 
JJ=J-1 
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1 A WRITE (6,16). JJ, {DUMMY I (J, I) , 1 = 1,6) 

16 FORMA I <2X, I3,F16.A,F15.A,F12.1,FL3.3,F1A.A,F15.A,F19.1, F16.2) 

WRITE 16,181 

18 FORMAT! 1H1///1X,9H SEC. NO. ,2X,9HGAS CQEF . , 2X, lOHCOOL .COEF. , 2X, 

1 10HAREA C0RR.,3X,9HR0UGHNESS,3X,llHRAD.0F CURV , AX, 6HQ RAD., 

25 X, 8HNUC • HEAT , AX , 12HSHELL* THICK »,2X, 16HCDAT .THK.OR RES.//50X, 5HC IN 
3. ), BXoHi IN. ),7X, 8HBTU.PER.,3X,BHBTU.PER. , 7X* 5H( IN. ) , 6X, 6H! IN) OR 
A, llHSu. IN-SEC-R/75X, 9HSQ. IN-SEC, 2X, 9HCU. IN-SEC,25X, 8HPER.BTU. /// ) 

00 20 J=2, NXX 
J J = J- 1 

20 WRITE! 6,22 ) J J , t DUMMY I ! J , I ) , I =9, 1 7 > 

22 F0RMAT(3X,I3,2X,2F11.A,F12.A,E15.3,E1A.A,F11.3,2F13.3,F17.A) 

WRITE! 6, 2A ) DEN , E S , AMU, BRA Z 

2A FORMAT I///1X, 18 HA VG.METAL OEMS ITY= , F6. 3, 2X, 9HLBS/CU. I N, 3X, 

1 22HSHELL YOUNG'S MODULUS 2 *, E10. 3, 2X»3 HP SI»3X,i6H POISSON'S RAT 

2IO=,F6.3,3X, 12HBRAZE TH ICK= , F7. A, 2 X , 5H! I N . 1 ) 

WRITE! 6,26) WG, WL , TGT, TL I N, PGS ! NXX ) , PL I N, D I AM, TLENG , ARCON , ARD I V 
26 FORMAT l 1H1//9X, 21HPR0PELLANT FLOW RATE , F9. A, 2X, 7HL0S/ SEC , 23X, 

1 17HC00L AN T FLOW RATE , AX, F9. A , 2X , 7HLBS/ SEC// 9X , 1 9HCHAMBER TEM 

2PERATURE, 3X,F8. 1,2X,6H0EG.R , 2AX, 22HENTRANCE COOLANT TEMP., 
3F8.2,2X, 6HDEG. R //9X, 16HCHAMBER PRE SSURE , 6X, F8 . 1 , 2 X , 6HPS I A , 
A2AX.22HENTRANCE COOLANT PRESS, F8 . 1 , 2X , 6HPSI A /// A3X , 15HI HRUAT DIA 
5METER, 13X,F9.A,2X,5H! IN. )/A3X, 19HTUTAL NOZZLE LENGTH, 8X, F 10. A , 2X, 
65H! IN. )/A3X,26HC0NVERGtNT AREA RATIO A/ A* , 3X, F8 . A , / A3X, 25HD I VERGEN 
7 T AREA RATIO A/A«,AX,F8.A//) 

IF! LESTC1 ) GO TO 30 
WRITE! 6,28 ) 

28 FORMAT ! A2X, 1 6H C-l IS NOT USED/) 

GO TO 3A 
30 WR I TE ( 6, 32 ) 

32 FORMAT l A2X , 1 6H C-l IS USED /) 

3A IF(LEsTC2l GO TO 38 
WRITE! 6,36) 

36 FORMAT! A2X, 16H C-2 IS MOT USED/) 

GO 10 A2 
38 WR I TE ! 6 , AO ) 

AO FORMAT! A2X, 16H C-2 IS USfcD /) 

A2 IF!L£iTC3) GO TO A6 
WRITE! 6,AA) 

AA FORMAT! A2X, 16H C-3 IS MOT USED/) 

GO TO 50 
A6 WRITE!6,A8) 

A8 FORMAT! A2X, 16H C-3 IS USED /) 

5C IFtLESTCA) GO TO 5A 
WRITE! 6,52) 

52 FORMAT ^A3X,21HTUBC IS ASSUMED ROUGH/) 

GO TO 58 
5A WR 1 TE ! 6, 56 ) 

56 FORMAT! A3X,22HTUBE IS ASSUMED SMUOTH/) 

58 1FILESTC5) GO TO 62 
WRITt! 6,60) 

60 FORMAT ( A3X • 3 AHAS SUMEO FLAT PLATE HEAT CUNOUCT ION/ // > 

GO TO 66 
62 WR 1 Tfc I 6 , 6A ) 

6A FORMA T { A3X, 30HA SSUME 0 RADIAL HEAT CONDUCTION///) 

66 STORE! 8)=0»5 
STORE! 1 5 ) =0 • 25 
STORE! 16 ) = A . 0 
68 ALSEC = 0.0 
RADI - 0.0 
HE 1 GHT =0 • 0 
HE1GH1 = 0.0 
PEST=PLIN 
T W1 =0. 0 
TGCER=0.0 

1 T= 1 

ATOT = 0 .0 
UTUT=0.0 
C TRT=0 . 0 
CTRH=0 .0 
C TRG=0 . 3 
C TRR=0 . C 
CTR I=G .0 
C TRQ=0 . 0 
KKK = 1 
KK= 1 
K = 1 

70 DO 72 L=1 , 1 8 
JMT+1 

72 DUMMY I ! 1 , L ) = DUMMY I ! J , L ) 

IF! IT.EQ.l) GO TO 7A 

AL^EC = SQRT ( ( (DHG-DHG1 )*0.5)**2 + !AX- AX 11**2) 

7A ANGBRD=PHI*A2 
COSA=COS(ANGBRO) 

IT=IT+1 

IF(TESTC.NE.O.O) GO TO 76 

T OVERK =TKCER 

TKCER=0.0 
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c 


CALCULATION OF GEOMETRY 


76 WEBT=2«0*TW+BRAZ 

GGAS*WG/UOHG*»2)*PI/4.0) 

I F ( IT. GT. NUMB) GO TO 80 
BNMEB-ANWEB! IT+1) 

IF( BNWEB.GE. ANWEB) GO TO 80 
PHI I=PI/BNWEB 

RADIUS=(PHII*(DHG*0.5+( TW+TKCER) *COSA)-WEBT*0. 5) / l 1.0-PHI l*COSA> 
HOLD=2 *0*AR£AL-P 1 2*RADI US**2+RADIUS*TW 
HEIGHT=HULO/ (2.0* RADI US -TW) 

78 RADSH= ( HE I GHT+RAD I US+T W ) *COSA+DHG*0 . 5 
WIDTH=2.0*PHIi*RAO$H-WEBT 
T HIGHT =HEIGHT 

HElGHT=H0LD/{RADiUS+W!DTH/2.0-TW) 

XX=AB$ ( HE I GHT-THIGHT } 

CTRG=CTRG+1.0 
IFICTR6.GT .10.0) GU TO 84 
IF ( XX . GT *CT( 2) ) GO TO 78 
HIGH= HEIGHT* RADIUS 

WP=P I2*RADIUS+2« 0*HE IGHT+WI DTH/2* 0-TW+RAD IUS 
RADIUS=RAOIUS/2.0 
GO TO 88 
80 PH I 1=P I /AN WEB 

RADlUS=(PHII*(DHG*0.5+( TW+TKCER ) *COSA ) -WEBT*0. 5 ) / ( I . 0-PHI l*COSA) 
HOLD* AREAL-PI 2* RADI US* *2 
HE I GHT =HOLO/(2.0*RADIUS) 

82 RADSH- ( HE I GHT+RAD I US+TW) *COSA+DHG*0 • 5 
HIDTH=2.0*PHll*RADSH-WEBT 
THIGHT=H6IGHT 

HE IGHT=HOLD/ (RADIUS* WIDTH* 0.5) 

XX=ABS( HE I GHT-THIGHT ) 

CTRG=CTRG+I . 0 
IF(CTRG.GT.IO.O) GO TO 84 
I F ( XX. GT »CT( 2 ) ) GO TO 82 
HIGH=HEIGHT+RAOIUS 
WP=PI*RAQlUS+2. 0* HEIGHT* WIDTH 
GO TO 88 
84 WRITE(6,86) 

86 FORMAT t///22X,36HITERATI0N LIMIT EXCEEDED ON GEOMETRY///) 

GO TO 268 
88 CTRG=0 .0 

RAD=RADIUS+0.5*TW 

DHL=4.0*AREAL/WP 

AWALL*PI2*ALSEC*ANWEB*(RADI+RADIUS+TW+TWI)*AREACR 

RRATSQ=UA0IUS/RC)**2 

IF( LESTC5) GO TO 90 

RADO=RAD IUS 

RAOT W=TW 

IF(TEsTC.LE.O.O) GO TO 92 
R ADC=RAU IUS 
RADCER-TKCER 
GO TO 92 

90 RADO=RADIUS+TW 

RADTW=RADQ*ALOG( RADO/RAOIUS) 

IF(TESTC.LE.O.O) GO TO 92 
RADC=RAOO+TKCER 
R AOCER=RADC*ALOG ( RADC/RAOO) 

92 AREAL=AREAL*ANWE8 
GL I Q=WL/AREAL 
WLARDH=GLlO*DHL 
T ST AT 4=TGS 
PST34=PGS 
CALL STATE! 4) 

RG0=A4/RSTF4 

IF( IT.NE.2) GO TO 94 

FPLIQ=0.0 

QNUC-0.0 

QNUCPT=0 . 0 

GLIQA=GLIQ 

GO TO 96 

94 A VGWP= ( WP+WP1 ) /2 . 0 

AREALM={ AREAL+AREAL1 ) /2.0 
DHLAVG=4.0*AREALM/ANWEB/ A VGWP 
GLIQA=WL/AREALM 

QNUC=QGEN*SHELTH*ALSEC*PI*{RADSH+RAOSHl+{ SHELTH+SHELT1 ) *0. 5) 
QNUCPT=ONUC/( ALSEC*PI*{ RADSH+RADSH1 ) - WEBT*ANW6B ) 
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CALC FOR 


TEMP JF 


THE 


GAS MALL 


1 

\ 


h 


96 1 F ( TEiTC) 1J0,98, 100 
98 TGCER=QEST* TUVERK+TWGE 
GU TO 108 

100 AKCEK=KMET(TCER,5) 

XTCER=TCER 

TGcER=(OEST/AKCER )*RAUCER+TWGE 

TCER=(TGCER+TWGE)/2.0 

XX=A8S( ICER-XTCER) 

1F( XX.LT.CTt 3) > GO TO 106 
CTRT=CfRT+1.0 
IFICTRT-25.0) 103,132,102 
102 WRITEC6.106) 

106 FORMAT I///22X»66HITERATI ON LIMIT EXCbEOEU ON AVG COATING TEMP///) 
GO TO 268 
106 CTRT=0.0 

108 TGE=C TGT-TGS)*AKR+TGS 
IF(TESTC)112, 110,110 

110 TGF=TGS+0.5*(TGCER-TGS)+0.22*< TGE-TGS) 

GO TO 116 

112 TGF=TGS*0.5*{TWGE-TGS) +0.22*1 TGE-TGS) 

116 T STAT6 = TGF 
PST36=PGS 
CALL STATE ( 6 > 

CPGF=CPSTA6 

AMUGF=AMUST6*A 

RGF*A6/RSTF6 

AKGF=AKSTA6*A 

PRGF= ( CPGF*AMUGF) /AKGF 

T AKR = AKR 

AKR=PRGF* *0.33333333 
XX=ABS(AKR-TAKR> 

CTRR = L TRR+ 1.0 

I F I C1RR.GT.10.0) GO TO 116 
1F(XX.GT.CT(5>) GO TO 108 
C TRR = 0 .0 
GO TO 120 
116 MR I Tb ( 6, 1 1 8 ) 

110 FORMAT { // /22X,63HITERAT ION LIMIT EXCEEDED UN RECOVERY FACTOR///! 
GO TO 268 

120 REGF= 1 ( GGAS*OHG ) / AMUGF ) * RGF/RGO 
REEXU=REGF**.8 
122 PRGEX=PRGF**.3 

HGF ~ ( (PRGEX*COEFG*AKGF)/DHG)*REEXG 
IFITE6TC) 128,126,126 
126 TWG=TGE-( QLST-QR AD 1 /HGF-U t ST*TUVERK 
GO TO 130 

126 TWG=TGG-(QEST-URAD)/HGF-QEST*RADCER/AKCER 
GO TO 130 

128 TMG=TGE-(OESr-QRAU!/HGF 
130 TERMX=A8 SI TMGE-ThG) 

ABS1S=TMGE 

IFITERMX.LT. CTt 8) ! GD TU 162 
IFIIWG-TL! 132,132,136 
132 U=UEST*0.1 
GO TO 226 

136 IFtCTRI .GT.O.O! GO TO 136 
TMGE=t TMG+TWGE) *3.5 
GO TU 138 

136 CALL 1TKAT(TMGE,THG,TVIGXI,TWGY1,THGE) 

138 T WGX 1 = A8S I i 
TWGYl=TWG 
CTRI=CTR I +1. 0 
IFICTR1.LT.25.0) GU TO 96 
MR I TE ( 6 , 160) 

160 FURMA T (///22X, 66H I TER AT I ON LIMIT EXCEEDED ON TEMP OF GAS WALL///) 
GO TO 268 
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c 


CALC. OF TEMP. OF LIQUID WALL 


142 CTR 1*0 .0 

IF!TESTC.GT.O.O) QEST«QEST*RADC/RAOO 
144 AKM-KMET ( TMET1, 1 I 

TWL*TWG-( IUEST/AKM)*RADTW) 

TMET=(TWG+TWL)/2.0 
IF! TMET-TL ) 146,146,148 
146 Q«UEST*Q.l 
GO TO 224 

148 TERMX*ABS!TMET1-TMET1 
CTRl*CTRl+i.O 

IFICTR1.GT.25.0) GO TO 150 
IF(TfcRMX.LT.CT(3) ) GO TO 154 
TM6T1*TMET 
GO TO 144 
150 WRITE! 6, 152 ) 

152 FORMAT! ///22X»44H ITERATION LIMIT EXCEEDED ON TEMP OF LIQ WALL///) 
GO TO 268 
154 CTft 1*0 .0 

QSEC* l ( QEST+QNM1 ) /2 .0 ) *AWALL+QN JC 
QSECWL*QS£C/WL 

C CALCULATION OF THE COOLANT 

IF( IT .NE.2) GO TO 162 
IF(KK.GT.l) GO TO 206 
KK=2 
C3=1.0 
PL = P 

DELPM*C.C 
DELPF*0.0 
CFL IQ*0.0 
PST34-PL 
TSTA3*TL 
CALL STATE! 3 ) 

RQS*A4/R$TA3 

ROST*RQS 

156 PLTQT*PL+GL IQ*GL IQ*ROST /!ROS*RQ 5*772.176) 

PST34=PLTQT 

TSTA3*TL 

CALL STATE! 3 J 

HT=ENTH33 

ROT*A4/RSTA3 

HS=HT-A3*( GLIQ/ROSI **2 

ENTH33*HS 

PST34*PL 

CALL STATEI-3J 

EROS*ROS 

R0S=A4/RSTA3 

ERUST=ROST 

ROST= ( R0S+R0TJ/2.0 

XX=ABS ( ROS-EROS ) 

XXX=AtJSlROST-EROST> 

YY*ROST*CT ! 9 l 
CTRI*CTKI+1.0 

IF(CTRI.GT.IO.O) GO TO 158 
IF(XX.GT.YY) GO TO 156 
IF! XXX. GT . YY ) GO TO 156 
CTR 1*0.0 
TLSTAT*T$TA3 
SONIC V*S TOR El 9) 

GO TO 206 
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158 WR I TE( 6, 1601 

160 F0RMATI///22X, 54HITERATI0N LIMIT EXCEEDED ON INITIAL COOLANT CONDI 
1 T IONS// / ) 

60 TO 268 
162 DELTH=gSECWL 

IF(CTRQ.EQ.O.O) HS=HS1+DELTH 
ENTH33=HS 
PST34=PL 
T STA3=TLSTAT 
CALL oTATEl-31 
ROS=A4/RSTA3 
HT=HT 1+DELTH 
164 H£ST=HS 

HS=HT-A3MGLIQ/R0S)»*2 

XXX=ABS!HEST-HS) 

YY=ABSIHS*CT(6) ) 

C TRH=C TRH+ 1 • 0 

IFICTRH.GT.15.0) GO TO 16B 
IF! XXX.LT. YY) GO TO 172. 

166 PST 34=PL 
ENTH33=HS 
CALL STATEI-3I 
R0S=A4/RSTA3 
SON ICV=ST0RE(9) 

TLSTAT=TSTA3 
GO TO 164 
168 WRITE I 6* 1 70) 

170 FORMAT I///22X f 4 3H ITERATION LIMIT EXCEEDED ON STATIC ENTHALPY///) 

GO TO 268 
172 CTRH=0.0 

DELPM=GL I QA* ( GL I Q/ROS-GOVE RR ) /386.088 

PST34=!PL+P>/2.0 

ENTH33= I HS+HS 11/2.0 

CALL STATE ( -3 ) 

TLAVGS=TSTA3 
R0AVG=A4/RSTA3 
TSTAT4=TSTA3 
CALL STATE! 4) 

AMULAV=AMUST4*A 
RESTAV=GL I UA*DHL AVG/AMULA V 
IFCLESTC3) GO TO 174 
C 3= 1 . 0 
GO TO 176 

174 C3=!REsTAV*RRATSQ)**.05 
IFIC3.LE. 1.0937) C3 * 1.0 
176 F2=F 1 

IF! Lb S TC4 ) GO TO 178 

FFLIg=! 1. 0/ ( -4. C * AL0G1C ( ROUGH/DHLA VG/3. 7+1. 255/REST A V/F1» *0.5) ) ) 
1**2 

GO TO 18G 

178 FFLIO=! 1*0/1 4, 9» AL0G1QI F1**0.5*RESTAV)-0.40) )**2.0 
180 F1=FFLI0 

C TRR = C TRR+ 1.0 
IF(CTRR.GT.IO.O) GO TO 182 
DELTF=ABS! F1-F2) 

IF! DELTF.GT.CT! 1) ) GO TO 176 

C TRR = 0 • 0 

CFLlg=FFLlO*C3 

OELPF=GLIQA*GLlOA/ROAV6*CFLIQ*ALSEC/DHLAVG/193.044 
P EST = PL 

PL=P-UtLPM-DELPF 
ABS I S=PL 

IF!PL-PLIN*C.3) 186, 186, 190 


85 



182 MR I TE! 6t 184) 

184 FORMA T(///22X t 43HITERAT ( DM LIMIT EXCEEDED ON FRICTION FACTOR///! 
GO TO 268 

186 PLIN=PLlN+50*0 
P = PL IN 
TL = T1 
QE S T = U I 
QNMl =Q I 
WRITE! 6, 188) 

188 FORMAT! ////22X,28H, NOZZLE COOLANT PRESSURE* LOST////! 

GO TO 68 

190 XX=ABS(PEST-PL) 

IF(CTR1 .GT.25.0) GO TO 192 
IF! XX.GT.CT(4) ) GO TO 196 
C TR I = 0 • 0 
GO TU 200 

192 WR ITE( 6f 1 94) 

194 FORMAT !///22X f 49HlTERAT ION LIMIT EXCEEDED ON L III PRESSURE///! 

GO TO 268 

196 IF(CTRI.EO.O.O) GO TO 198 

CALL ITRAT! PEST # PL , P 3X , PQ Y, PL I 

198 PQX=PEST 

pgy=ABSis 
CTR I=CTR I + 1 • 0 
GO TO 166 

200 PST34=lPLT0Tl+PLT0T)/2.0 
ENTH33=(HT+HT1 ) /2.0 
CALL STATt(-3) 

R0AVGT=A4/RSTA3 

TLAVGT=TSTA3 

DELTPr = ROAVGT*Al»U£LTH*( TL A VGT-TLA VGS ) /TL A VGS*DELPF*ROAVGT 
1*TLAVGT/TLAVGS/R0AVG 
PTfcST=PLTUT 
PLTOT = PLTOT l'-OELTPT 
XX=ABS!PTEST-PLTOT) 

CTRI=CTRI+1.0 

IFICTRI-GT.15.0) GO TO 202 

IF! XX.GT.CK4) ! GO TO 200 

CTR 1=0.0 

P $T34=PL TOT 

ENTH33=HT 

CALL STATE! -3) 

TL=T3TA3 
R0T=A4/RSTA3 
GO TO 206 

202 WR I TE ( 6 1 2 04 ) 

204 F0RMAT(///22X,46HITERATIQN LIMIT E XCEEOEL) ON LIQ TOTAL PRESSURE//) 
GO TO 268 

206 TSTAT4=TLSTAT 
PST 34=PL 
CALL STATEl 4 ) 

AMUBLK=AMUST4*A 

R0S=A4/RSTF4 

RE8ULK = WLARL)H/AMU8LK 



C CALC. UF HEAT FLUX 

IF(THL-TL) 208,238,210 
208 Q=«JEST*0.1 
GO TO 224 

210 TLF = { TL + TWO/2.0 
T ST AT 4=TWL 
CALL STATE (4) 

AMUTtfL=AMUST4*A 
RUTWL=A4/RSTF4 
T ST AT4 = TLF 
CALL STATE ( 4) 

CPLF=tPSTA4 
AHULF=AMUST4*A 
AKLF = AK.STA4*A 

RELF=WLAR0H*A4/IRSTF4*AMULF*R0S) 

IF(L£STC2) GO TO 212 

C2=1.0 

GO TU 214 

212 C2=(RfcLF*FtRATS0)**.05 
IF(C2. Lb. 1.0937) C2=1.0 
214 I F ( LES TCI ) GU TO 216 
Cl=1.0 
GO TO 218 

216 ZOVZ= l AMUT tfL/ROT WL ) / ( AMU8LK/R0S ) 

Cl=1.0+( . 01457* ZD VZ ) 

218 PRLF=(CPLF*AMULF)/AKLF 
REfcXL=RfcLF**.8 
PRbX4=PRLF**.4 

P AHLF= { PREX4*C0EFL» AKLF /UHL ) *REEXL 
1F(RC 1220*220,222 
220 C2=1.0 
222 HLF=PAHLF*C2*C1 
Ul={TrfL-TL)*HLF 
Q=Ul*RAUl US/ R ADO 
224 AGS 1 S = uEST 

T ERMX= ABS (OEST-Q) 

IF[TfcRMX.GT.CT( 7) ) GO TO 226 
C TR0 = 0 • 0 
GO TO 236 

22 6 IF(CTRU.GT.O.O) GO TO 228 
gEsr=(U*QEST)*0.5 
GU TU 230 

22B CALL 1TRATJ0EST, J, TOX 1 , TUY 1 , JEST ) 

230 T UX 1 =AdS I S 
T UY 1 = 1/ 

IFITtSTC.GT.O.O) JEST=UEST »RADO/R ADC 
CTRI=0.0 
CTRQ=CTRU+ 1.0 
I F ( C f KQ-2 5 . 0 ) 96,232,232 
232 WRITfc{6,234> 

234 FORMAT t///22X,37HITERATI UN LIMIT EXCEEDED UN HEAT FLUX///) 
GO TU 26B 

236 VEL=( GL lu/ROS) *A 
AMACH = VfcL / SON I C V 



C CALCULATION OF STRESS 

TMETA=tTMET+TL)/2#0 
ALPHA1=KMET ( TMETA , 4) 

ALPHA2=KMET(TMET, A) 

E-KMETtTMET, 3) 

S IGALM=KMET< TMET, 2 ) 

IFl IT.NE.2) GO TQ 238 
SNS = 0.0 
GO TO 240 

238 SNS=-0.125*(PGSi+PGS)*l0HGl**2-DHG**2l/(UHG*C0SA) 
1 +SnS*DHG1/DHG*C0S(ANG8RD-ANGRD1 ) 

240 SNPSI=PGS*DHG*0.5/COSA 

EPSLNM=(SNS-AMU*SNPSI >/<ES*SHELTH) 

EPSLNH=( SNPS I-AMU*SNS ) / ( ES*SHELTH) 

TERM1MPL-PGS)*RAD/(E*TM) 

EPSLT=-TfcRMl*AMU+ALPHAI*( TMET-TL) 

DELT=( ( TERMl+ALPHAlM TMET-TL) l-AMU* ( EPSLNH-EPSLT ) )*2.0*RAD 
i-(PI*EPSLNH/ANMEB)*(DHG«-2.0*RAD/C0SA) 

TERM2=fc*ALPHA2*( TWG-T ML ) / C 2 . 0* { 1.0-AMU) ) 

TERM3=i.0681*E*TW*DELT/( { 1 . 0-AMU**2 ) *RAD**2 ) -TERM2 
TERM4=TERMl*E-0.2796*E*TM**2*DELT/( <1.0-AMU**2)*RAU**3) 

S IGGUT=TERM4+TERM3 

SIGIN=TERM4-TERM3 

EPSLTX=EPSLNM-EPSLT 

SlGLON=E*EPSLTX 

SIGL0I = S1 GLDN + TE RM2 

SIGLUU=SIGL0N-TERM2 

ANGR01=ANGBRD 

PGS1=PGS 

IFl IT.NE.2) GO TO 242 

WGHT=0.0 

GO TO- 244 

242 SHELTA=I SHELT1+SHELTH 1/2.0 

MGHT=WGHT+( PI2*ALSfcC*SHELTA*DEN)*(DHGH-DHG+2iO*<HEIGHl+ 
IHEIGHT+RAD1+RADIJS+TW+TM1+SHELTA)*CQSA)+{ALSEC*ANMEB 
2*0. 5* ( I M+TM1 ) *OEN ) *(PI2*(RADl + RADIUS)+rlEI GH1+HE I GHT ) 

C OUTPUT AND NEXT STATION 

244 P=PL 

ROVT=RAD/TW 
ATOT=ATOT+AWALL 
QTOT=QTOT+QSEC 
I F ( K . EQ . 2 ) GO TO 250 
IF(TESTC) 248,248,246 
246 IFITCER.GT.TEMP3) K = 2 
248 IF l TMET .GT .TEMP 2) K=2 
IF(TMETA.LT.TEMPI) K=2 
250 IFl IT. GT. NUMB) GO TO 258 
HTI=HT 
HSl=HS 

GOVEKK=GL IQ/ROS 
PLIOT 1=PLT0T 
0EST=Q 
QNM1-U- 
DHG1=0HG 
A XI =AX 

HE IGHi=HElGHT 
TM1=TW 

SHELT 1=SHELTH 
RAOSHi=RAOSH 
AREAL I=AREAL 

IFl BNMtB-ANWEB) 252,256,254 
252 WP1=2.0*( MP-HEIGHT-RADIUS+TW) 

RAD1=2.0*RADIUS 

K.KK=2 

GO TO 258 

254 WPl=WP/2.0+HEIGHT+RADIUS-TW 
RADI = RAD I US/ 2.0 
KKK = 2 
GO TO 258 
256 MP1=MP 

RA01=RAJ1US 

258 IF ( IT.NE.2) GO TO 262 
MR 1 TE ( 6 , 260 ) 

260 F UR MATI2X, 7HSEC.N0. ,2X,7HQ/A 0UT,4X,6HU/A IN , 3X, 9HEFF . GAS T,2X, 

1 6HT. REF., 3X,8HT. COAT. G,3X,7HT .MET .&, 3X, 7HT. MET. L, 5X,6HL IO.T. 

2, 4X,6HLIM.P. ,4X,7HLIU.VEL,4X,7HMACH N0/84X, 5HT0TAL , 5X , 6HSTAT IC/ 
310X,18rt(BTU./SQ. I N .-SEC . ) , 4X, 7H( DEO . R ) , 3X , 7H ( DEG. R ) , 4X, 7H ( OtG. R ) , 
43X,7H(UEG.R),3X,7H(DEG.R),4X,7HIDEG.R) ,3X,6H(PSIA) , 4X.6HFT/SEC/// ) 
262 DO 264 L*I,39 

264 DUMMYUt IT,L) =DUMM YO ( 1 ,L ) 

J J 3 I T- 1 

MR I T t ( 6 , 266 ) JJ,Q,UI,TGE,TGF,TGCER.TMG,TWL,TL,P,VEL,AMACH 
2 66 FORMAT (3X, 13, 2F 1 1 . 2 , 5F10.0, FI 1 . 1 , 2F 10. 0 , F 12 . 3 ) 

I F( NUMB.GE .IT) GO TO 70 
GO TO 270 
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FINAL OUTPUT DATA 


t 


c 

c 

c 

c 

c 


268 NXX= I T**l 

270 IFUKK.EO.l) GO TO 274 
WRITE (6,272) 

272 FORMAT !///25X,4LH****N0TE NOZZLE DESIGN CONTAINS A SPLICE//) 

274 IF(K.fcU.l) GO TO 278 
WRI TE ( 6, 276) 

276 FORMAT (///25X»60H***» NOTE TEMPERATURE RANGE EXCEEDED ON MATERIAL 
I PROPERTIES) 

278 WR i TE ( 6, 280) 

280 FORMAT ( IHI///2X, 7HSEC . NO. , 2X, IOHFR I .P . DROP, 2X , 1 OHMOM. P. DROP , 2X, 
i 6HFR I ■ FACT , 7X » 3HHGF , 1 1 X, 3HHLF , IOX , 2HC 1 , 9X , 2HC2 , 8X , 2HC3, 7X , 

29HU/SECT ION, 6X, 7 HQ/ A NUC// 1 3X, 5H< Pi I ) , 8X , 5H( PS I) , 1 7X, 2 lH( BTU/SO . I N 
3. SEC. DEG. R) ,39X, 7HBTU/SEC , 4X, l 3HBTU/ SQ. I N . SEC/// ) 

DO 282 J=2» NXX 
J J = J-i 

282 WR I TE ( 6, 284 ) J J , ( OUMMYO! J , I I , I *1 , ID ) 

2 84 F0RMAT(3X, 1 3 , 2F 1 2 . 2 , 3F 1 3. 5 , 3F 1 1 . 3 , F 12 . 1 , F 1 5 . 3 ) 

HR ITE ( 6, 286) 

286 FQRMAT(irtl///2X, 7 H SEC. NO. , 6X , 9HSHELL RAD, 8X, 

1 11HTUBE RADIUS, 6X,IIHTUBE HE IGHT »6X» 10HTUBE WI DTH, 9X, 7HHY0. D 

2 1 A, 7X , 9HWALL ARE A , 6X , 1 7HT0TAL TUBE HE I GHT// 18X , SHI I N. ) , 1 i X, 5H! I N. ) 
3, 12X.5HI IN.), 13X,5H( IN. ) , 1 2 X, 5H( I N . ) , 1 DX , 7H( SQ . I N ) , I2X, 5H ( IN. )///) 
DU 288 J =2 , NXX 
J J= J-l 

288 HRITE( 6, 290) J J , 1 OUMMYOt J, I ) , I =11 , 17) 

290 FORMAT (3X, I3,7FI7.4) 

WRITE! 6,292) 

2 92 FORMA I ( 1 HI ///2X , 7HSEC . NO. , 2 X, 6HL I Q . T. , 3X , 6HL I Q. P. , 6X , 

1 I2HLIQ. ENTHALPY, 1IX, 11HLIQ. DENSITY, 9X.6HRE GAS,5X,6HPR GAS, 

24X , 6HRE L 1 U, 7X, 6HRE LIg,6X,6HPR L I Q/U X , 6HSTAT I C , 3 X , 5HT0T AL , 4X, 6HS 
3TATIC,6X,5HT0TAL,6X,6HSTATIC, 5X, 5HT0TAL , 7X , 4HF I L M , 7X, 4HF I LM, 6 X, 
44HBULK,9X,4HFILM,BX,4HFILM/11X, 7H(DEG.R) , 2X , 4HPS I A, 5X , 6HBTU/LB, 6X , 
56H8TU/LB, 9X, 9HLBS/CU. IN///) 

DO 294 J=2 , NXX 
JJ=J-1 

294 WR I TE ( 6,296) J J , ( DUMMYO ( J , I ) , I = l B , 28 ) 

296 FORMAT ( 3X, 1 3 , F 10 . 1 , F9. 0 , 2F 1 1 . 2 , 2f I 1 .5, E 13 . 3, F9 . 3, 2E 1 3 . 3 , F9 . 3 ) 

HR 1 TE ( 6, 298 ) 

298 FORMAT! 1HI///IX, 7HSEC.ND. , IX, 8HR OVER T,2X,I0HTUBb DELTA, 22X, 

I I8HTUBE WALL STRb SSE S , 1 5X, 9HA VG . ME TAL , 3X , SHY I ELD, 5X , 

2 7 HCUEF . UF , 6X, 7H YOUNG ' S/ 31 X , 7HL0NG. A V, 3X, 8HL0NG . OUT , 4X , 7HLUNG. IN, 
34X,8HTANG.QUT,3X,7HTANG. I N , 4X , 5HTEMP. , 3X, 8HSTRENG TH, 4X, 7HLIN.EXP, 
46X»7HM0UULUS/21X, 5H{ I N . ) , 6X , 5H! PS I ) , 6X, 5H ( P S I > , 6X, SHI PS1 ) , 6X, SHI PS 
51) ,6X,5H! PS I ),5X, 5H0EG.R,4X,5HIPSI) f 6X,7H( IN/ IN ) , 7X, 5H ( PS I ) // / 1 
DO 3G0 J=2 , NXX 
JJ=J-1 

300 WRI TE ( 6, 302 ) J J , ( DUMMYO! J , I ) , 1 =29 , 39) 

302 FORMAT! 2 X, 1 3 , FI 0 . 3 , F 11 . 5 , 5F 1 1 .0 , F9 . 0» FI 1 . 0, ZE 1 3 . 3 ) 

DELTT=TL-TLIN 

DELTP=PLIN-PL 

WRITE (6,304) g TOT , ATOT , HGHT , DELTT , UELTP 
304 FORMAT! 1H1////9X, 27HTOTAL COOLANT HEAT PICKUP = , F 1 0 . 2 , 2X , 7HBTU/SEC 
1 / / 9X » 27HT0T AL NOZZLE WALL AREA = , F 10 . 2 , 2 X , 6HSQ. IN.//9X, 

227HTDTAL NOZZLE WEIGHT = , F 1 0 . 2, 2X , 3HLBS// 9X , 2 7HT0TAL CUOLANI 

3 TEMP . DtL T A = , F 1 0 . 2 , 2X , 5HDEG. R//9X , 27HT0T AL COOLANT PRESS. DELTA = 
4,F10.2,2X, 3 HP SI ) 

WRITE! b, 306) 

306 FORMAT ( //SOX, 5HT I TLE//50X, 7HRESULTS/////50X, I6HFJTURE REVISIONS//) 
GO TU 6 
END 

SUBROUTINE 1TRAT( V, W, X, Y, Z ) 

IF! V-X)2, 1,2 

1 X=X*I.l 

2 OENOM=V-X 
SLUPE=! W-Y)/DENOM 
IFtSLUPE.EU.l.Q) GO TU 3 
UNMI-1. -SLOPE 

Z=ABS ( (X*SLOPE-Y)/DNMI ) 

GO TO 4 

3 Z*W 

4 RETURN 
ENO 

REAL FUNCTION KMETITX.N) 

COMMON X(10,5),TEMP1,TEMP2,TEMP3 
DIMENSION TEN! 5 ) 

DATA TEN/l.DE-4, 1.0E3,1.0b6, l.OE-6, 1.0E-4/ 

N = 1 » THERMAL CONDUCTIVITY 
N = 2, YIELD STRENGTH 
N=3, YOUNGS MODULUS 
N=4, THERMAL EXPANSION 

N=5, THERMAL CONDUCTIVITY OF CUATING 
T = TX 

IF(TtMPl.GT.T) T=TEMP1 
IF(N.EU.S) GO TO I 
IF! TEMP2.LT . T ) T=TEMP2 
GO TU 2 

1 I F ( TEMP3 . L T . T > T = TEMP3 

2 W/1000.0 

KMET=(X(l,N)+T»!X(2,N)+T»(X(3,N)+T«(X!4,N)+T»(X(5,N)+T+!X(6,N) 
1+T*(X17,NJ+T*< XI B,N)+T*(X(9,N)+T*XI 10, N) )))))) ) I )*TEN(N) 

RETURN 

END 
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